I 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


s4<y%3 


AN 

INVESTIGATION  INTO  THE  USE  OF  AN 
EXISTING  SHOCK  TUBE  AS  A  DRIVER 
FOR  A  HYPERSONIC  SHOCK  TUNNEL 

by 

Mi 

chael 

H.  Sherman 

March  1989 

Th 

esis 

Advisor : 

Raymond  P.  Shreeve 

Approved  for  public  release;  distribution  is  unlimited 


T242349 


REPORT  DOCUMENTATION  PAGE 


REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


lb    RESTRICTIVE  MARKINGS 


SECURITY  CLASSIFICATION  AUTHORITY 


DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


3     DISTRIBUTION /AVAILABILITY  OF   REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5    MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


NAME  OF  PERFORMING  ORGANIZATION 

aval    Postgraduate    School 


6b    OFFICE  SYMBOL 
(It  applicable) 

Code    67 


7a    NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


ADDRESS  (City,  State,  and  ZIP  Code) 
Hiterey,    California      93943-5000 


7b    ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Monterey,  California   93943-5000 


NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b   OFFICE  SYMBOL 
(If  applicable) 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


ADDRESS  (City,  State,  and  ZIP  Code) 


10    SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 
NO 


TASK 
NO 


WORK  UNIT 
ACCESSION  NO 


TITLE  (Include  Security  Classification) 
I    INVESTIGATION    INTO    THE    USE    OF    AN    EXISTING    SHOCK    TUBE    AS    A    DRIVER    FOR    A 

'PERSON IC   SHOCK  TUNNEL 


PERSONAL  AUT.HOR(S) 

lerman,    Michael    H. 


a    TYPE  OF  REPORT      . 

ister's    Thesis 


3b    TIME  COVERED 
FROM  TO 


14    DATE  OF  REPORT    (Year,  Month,  Day) 

1989^    March 


5    PAGE  COUNT 

152 


i    SUPPLEMENTARY  NOTATION 

]ie  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
Dlicy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


COSATI  CODES 


FIELD 


GROUP 


SUB  GROUP 


18    SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Shock  Tube;  Tailoring;  Shock  Tunnel 


ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by 

Experiments  were  carried  out  u 
ibe  connected  to  a  two-dimensional 
on  of  steady  reflected  pressure  f 
irough  tailored  operation  for  inci 
lire  Helium  and  a  70%  Helium/30%  Ni 
pspectfully . 

Spark  and  continuous  light  sha 
vtical  window  at  the  Mach  4 . 3  loca 
iiration  flow  phenomena  in  a  shock 
listing  equipment.   Calculations  s 
id  15  m3  (6 '  diameter  x  15'  long) 
•nic  facility  for  instruction  and 


block  number) 

sing  an  existing  tube  alone,  and  with  the 

wedge  nozzle.   The  range  of  maximum  dura- 
rom  3.5  to  5  milliseconds  was  achieved 
dent  shock  strengths  of  3.4  and  2.0,usina 
trogen  mixture  as  the  driver  gas 

dowgraph  techniques  were  attempted  using  an 
tion.   Results  demonstrated  that  the  short 
tunnel  can  be  recorded  successfully  using 
howed  that  the  addition  of  a  Mach  10  nozzle 
dump  chamber  would  provide  a  useful  hyper- 
research  . 


:    DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
0  UNCLASSIFIED/UNLIMITED       □  SAME  AS  RPT  □  DTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


a    NAME  OF  RESPONSIBLE  INDIVIDUAL 

rof.  Raymond  P.  Shreeve 


nbCibW^c-^riaCode) 


t.2c    OFFICE  SYMBOL 

Code    67Sf 


)FORM  1473,  84  mar 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Ou  s    Govtrrm  int  Printing  OMici     !»(«— COS  24. 

UNCLASSIFIED 


Approved  for  public  release;  distribution  is  unlimited 

An  Investigation  into  the  Use  of  an  Existing 
Shock  Tube  as  a  Driver  for  a  Hypersonic  Shock  Tunnel 

by 

Michael  H.  Sherman 

Civilian,  Naval  Air  Systems  Command 

B.S.,  Pennsylvania  State  University,  1983 


Submitted  in  partial  fulfillment  of  the 
requirements  for  a  degree  of 


MASTER  OF  SCIENCE  IN  AERONAUTICAL  ENGINEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1989 


ABSTRACT 

Experiments  were  carried  out  using  an  existing  tube 
alone,  and  with  the  tube  connected  to  a  two-dimensional 
wedge  nozzle.   The  range  of  maximum  duration  of  steady 
reflected  pressure  from  3.5  to  5  milliseconds  was  achieved 
through  tailored  operation  for  incident  shock  strengths  of 
3.4  and  2.0,  using  pure  Helium  and  a  70%  Helium/30%  Nitrogen 
mixture  as  the  driver  gas  respectfully. 

Spark  and  continuous  light  shadowgraph  techniques  were 
attempted  using  an  optical  window  at  the  Mach  4.3  location. 
Results  demonstrated  that  the  short  duration  flow  phenomena 
in  a  shock  tunnel  can  be  recorded  successfully  using 
existing  equipment.   Calculations  showed  that  the  addition 
of  a  Mach  10  nozzle  and  15  m3  (6*  diameter  x  15'  long)  dump 
chamber  would  provide  a  useful  hypersonic  facility  for 
instruction  and  research. 
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I.   INTRODUCTION 

The  simulation  of  flight  conditions  using  ground  based 
facilities  has  been  an  essential  part  of  all  aircraft 
development.   Figure  1  illustrates  the  various  types  of 
facilities  which  have  evolved  to  simulate  flight  at 
progressively  higher  speeds  and  altitudes.   Increasing 
energy  requirements  led  to  blow-down  wind  tunnels.   The  need 
to  avoid  condensation  in  nozzle  expansions  required  heaters. 
The  very  high  stagnation  enthalpies  associated  with 
hypersonic  flight  led  to  very  short  duration,  intermittent 
techniques.   Each  of  the  types  shown  in  Figure  1  was  in 
operation  more  than  20  years  ago.   A  renewed  interest  in 
hypersonic  simulation  has  been  generated  by  the  National 
Aerospace  Plane (NASP)  program. 

A  recent  review  of  facility  requirements  for  hypersonic 
propulsion  system  testing  by  Dunn  and  Lori  [Ref.  1] 
concludes  that  a  single  ground  based  facility  will  not  be 
capable  of  providing  the  necessary  testing  to  cover  the  NASP 
flight  envelope.   The  planned  use  of  Computational  Fluid 
Dynamic  (CFD)  codes  to  bridge  the  gap  left  by  inadequate 
ground  based  facilities  will  make  CFD  codes  an  essential 
part  of  the  design  process  but,  before  confidence  can  be 
placed  in  these  codes,  fundamental  experiments  are  needed  to 
validate  the  models  incorporated  for  real  gas  behavior. 
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Figure    1.      Ground   Based   Testing   Facilities   [Ref.    1] 


Recognizing  the  need  for  CFD  validation,  the  feasibility 
of  establishing  a  shock  tunnel  capability  by  modifying  the 
existing  Naval  Postgraduate  School  (NPS)  three  inch  diameter 
shock  tube,  was  examined  in  the  present  study.   The  time 
duration  and  magnitude  of  constant  pressure  available  in  the 
reflected  region  of  the  shock  tube  were  evaluated.   The 
concept  of  tailoring  was  examined  using  two  types  of  driver 
gases,  pure  Helium  (He)  and  a  mixture  of  70%  Helium  and  30% 
Nitrogen  (N2)  ,  driven  into  Nitrogen  as  a  test  gas.   An 
exploratory  shock-tunnel  nozzle  experiment  was  conducted, 
and  the  capability  of  recording  short-duration  flow 
phenomena  using  available  visualization  techniques  was 
demonstrated. 

It  was  concluded  that  the  development  of  an  operational 
hypersonic  shock  tunnel  was  feasible  and  would  provide  the 
means  for  follow-on  efforts  in  the  area  of  CFD  code 
validation.   At  the  same  time,  such  a  facility  would  provide 
a  useful  illustration  of  hypersonic  flow  to  aeronautical  and 
space  applications  engineering  students. 

The  following  material  is  divided  into  three  sections; 
1)  a  review  of  shock  tube  and  shock  tunnel  ideas,  2)  a 
discussion  of  the  experimental  investigation  including  a 
description  of  the  facilities,  experimental  procedure  and 
results;  and  (3)  an  evaluation  of  nozzle  throat  and  dump 
chamber  sizing  requirements  for  a  hypersonic  tunnel. 
Instrumentation  modifcations  will  be  required  in  the 


development  of  a  hypersonic  capability.   Appendix  A  reviews 
the  existing  facility  and  provides  recommendations  for 
upgrading  the  current  hardware  and  instrumentation. 
Simplified  flow  models  are  developed  in  Appendices  B  and  C 
to  enable  the  selection  of  nozzle  throat  diameter  and  dump 
chamber  volume  respectively.   Experimental  pressure  and  thin 
film  traces  obtained  in  the  experimental  program  are 
contained  in  Appendix  D  together  with  predicted  results 
obtained  using  a  computer  performance  code,  which  is  listed 
in  Appendix  E. 


II.   THEORETICAL  CONSIDERATIONS — A  REVIEW 

A.   BASIC  SHOCK  TUBE 

A  shock  tube  is  an  economical  short  duration  test 
facility  used  in  laboratory  experiments  in  gas  dynamics, 
aerodynamics  and  physics.   The  principles  of  shock  tube 
operation  ,  which  are  well  documented  in  numerous  texts,  for 
example,  Gaydon  and  Hurle  [Ref.  2]  and  Owczarek  [Ref.  3], 
are  illustrated  in  Figures  2  and  3  [Ref.  4].   As  shown  in 
Figure  2,  the  device  consists  of  a  long  tube  divided  into  a 
high  pressure  and  a  low  pressure  section  by  a  diaphragm. 
When  the  diaphragm  separating  the  high  pressure  driver 
section,  region  4,  and  the  low  pressure  section,  region  1, 
is  ruptured,  a  series  of  compression  waves  coalesce  into  a 
single  shock  front  (S)  which  compresses  and  heats  the  region 
1  test  gas  to  region  2  conditions.   Simultaneously,  a  series 
of  rarefaction  (expansion)  waves  (R)  travel  into  the  high 
pressure  section,  region  4,  cooling  and  expanding  the  driver 
gas  to  region  3  conditions. 

Region  3  (the  expanded  driver  gas)  and  region  2  (the 
compressed  test  gas)  are  separated  by  a  contact  surface  (C) . 
The  gases  on  either  side  of  the  contact  surface  travel  with 
the  same  velocity,  and  are  at  the  same  pressure,  but  are 
different  in  temperature.   A  time  history  of  the  process 
(time  vs.  distance)  is  shown  in  Figure  3. 
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Figure  2.   Thermodynamic  States  in  the  Operation  of 
a  Shock  Tube  [Ref.  4] 
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Figure  3.   Time  History  of  Events  after  Diaphragm 
is  Broken  [Ref.  4] 


What  is  shown  in  Figure  3  can  be  determined  accurately 
using  the  method  of  characteristics  and  shock  theory  (see 
[Ref.  3]).   As  the  shock  front  reflects  off  the  driven  tube 
endwall,  the  region  2  gas  is  effectively  shocked  a  second 
time,  creating  in  region  5  stagnant  gas  at  high  pressure  and 
temperature.   This  region  may  be  used  for  high  temperature 
gas  dynamic  studies  or,  when  coupled  to  a  nozzle,  made  to 
serve  as  a  reservoir  of  compressed  and  heated  gas  for  a 
hypersonic  wind  tunnel. 

B .   PERFORMANCE 

The  strength  of  the  generated  shock  (Ms)  is  related  to 
the  pressure  ratio  (PA/PX)  across  the  diaphragm  and  the  speed 
of  sound  in  the  undisturbed  gases  through  the  eguation  given 
by  Gaydon  [Ref.  2]. 
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Figures  4  [Ref.  5]  and  5  [Ref.  6],  which  illustrate  Equation 
(1) ,  show  that  the  strongest  shocks  are  generated  by  using 
high  pressure  ratios  and  by  using  driver  gases  with  high 
speeds  of  sound  and  low  ratios  of  specific  heats. 

As  the  desired  shock  strength  increases,  the  required 
PyPi  increases  significantly.  To  achieve  stronger  shocks 
several  approaches  are  possible.   The  most  common  include 
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Figure  4.   Pressure  Ratio  (P«/Px)  Required  to 

Generate  a  Given  Shock  Strength  Using 
Various  Gas  Mixtures  [Ref.  5] 
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(a)   Driver  Gas  Specific  Heat  of  1.4  and  a  Driven 
Gas  Specific  Heat  of  1.67 


Figure  5.   Variation  of  Shock  Strength  as  a  Function 
of  Speed  of  Sound 
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(b)   Driver  Gas  Specific  Heat  of  1.67  and  a  Driven 
Gas  Specific  Heat  of  1.67 


Figure  5  (CONTINUED) 
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the  use  of  an  area  reduction  between  the  driver  and  driven 
sections  [Ref.  5],  the  use  of  multiple  diaphragms  with  a 
buffer  gas  in  the  driver  [Ref.  5] ,  and  heating  the  driver 
gas  through  combustion  or  electrical  discharge  technigues 
[Ref.  6].   Figures  6  and  7  [Ref.  7]  show  the  relationship 
between  shock  strength  (Ms)  and  reflected  region  5 
conditions  for  eguilibrium  air. 

C.   SHOCK  TUNNELS 

The  reflected  region  5  pressure  and  temperature  may  be 
used  to  supply  a  supersonic  or  hypersonic  nozzle  if  the 
nozzle  throat  is  made  sufficiently  small  to  ensure  full 
reflection  of  the  incident  shock.   The  evolution  of  this 
arrangement  is  shown  schematically  in  Figure  8.   To  maximize 
test  time  and  minimize  the  nozzle  starting  time,  often  a 
diaphragm  separates  the  nozzle  and  dump  chamber  from  the 
driven  end  of  the  shock  tube.   The  dump  chamber  is  sized  and 
evacuated  to  a  level  which  ensures  that,  when  the  tube  is 
fired,  the  pressure  in  the  chamber  remains  below  the  static 
pressure  at  the  nozzle  exit  for  the  duration  of  steady  flow. 
[Ref  8:pp.  624-640] 

For  hypersonic  applications  a  conical  nozzle  is  often 
used.   Table  1  provides  calculated  test  conditions  for 
various  flow  expansions  using  ideal  gas  theory.   The 
computer  program  used  for  the  calculations  is  described  in 
Appendix  E.   Since,  in  an  isentropic  expansion,  static 
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Figure  6.   Reflected  Pressure  as  a  Function  of  Shock 
Strength  for  Equilibrium  Air  [Ref.  7] 
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Figure  7.   Reflected  Temperature  as  a  Function  of  Shock 
Strength  for  Equilibrium  Air  [Ref.  7] 
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TABLE  1 

THEORETICAL  PREDICTION  OF  TEST  CONDITIONS  FOR  A  HYPERSONIC 
NOZZLE  (Me  =  10)  DRIVEN  FROM  THE  REFLECTED  REGION 

OF  A  SHOCK  TUBE 


Driver 
Gas 

Ms 

P5 
(Atm) 

T5 
(°K) 

Pe 

(mm  Hg) 

(°K) 

~  2.4 

58 

973 

1.02 

50 

3.7 

32 

2000 

.55 

94 

Helium 

4.2 

23 

2460 

.41 

116 

4.5 

19 

2700 

.33 

131 

"l.8 

23 

632 

.4 

30 

Nitrogen 

2.4 
2.6 

9.7 
6.7 

971 
1110 

.169 
.118 

50 
52 

2.8 



5.3 

1200 

.09 

56 

Wtart 

(msec) 
.041 
.029 
.026 
.024 

.051 
.041 
.038 
.037 


temperature  decreases  as  the  Mach  number  (of  the  test 
section)  increases,  liquefaction  of  the  test  gas  becomes  an 
issue.   Since  nitrogen  closely  matches  the  aerodynamic 
properties  of  air,  but  liquefies  at  a  lower  temperature,  it 
is  preferred  for  hypersonic  facilities  in  which  simulated 
real  gas  effects  are  not  required.   This  is  illustrated  in 
Figure  9  which  also  shows  the  interesting  results  of  wind 
tunnel  experience  [Ref.  9].   In  practice,  condensation  is 
observed  to  occur  at  temperatures  which  are  lower  than  the 
saturation  temperature,  thereby  reducing  somewhat  the 
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stagnation  temperature  requirements  to  reach  a  given  Mach 
number. 

Given  the  extremely  short  test  time  available  in  a  shock 
tube,  any  coupling  to  a  hypersonic  nozzle  must  be  done  in 
such  a  way  so  as  to  minimize  nozzle  starting  transients  and 
provide  the  maximum  available  steady  flow  in  the  test  sec- 
tion.  The  starting  time,  defined  as  the  time  to  approach 
steady  flow  through  a  hypersonic  conical  nozzle  (7:  =  7/5)  , 
is  given  by  Smith  [Ref.  8]  as 


.      _  l!%  1/2  (1  -  ri)  +  4<Mp3/2-l)  +  oi(Me5/2-l)  ]  ( 2 , 

start   6  ac  L   e       M_    15  e        25  e 
5  e 


A  reduction  in  starting  time  may  be  achieved  by  using  a 
large  nozzle  angle  (3  =  nozzle  half  angle) ,  a  small  throat 
radius  (Rt)  and  a  high  stagnation  temperature  in  the 
reflected  region.   [Ref.  8:pp.  625-640]   The  variation  in 
starting  time  with  increasing  incident  shock  strength  for  a 
conical  nozzle  with  B  =  30  degrees  is  shown  in  Table  1. 

After  the  unsteady  starting  process,  the  useful  test 
flow  will  continue  as  long  as  a  sufficiently  high  pressure 
ratio  is  maintained  across  the  nozzle.   The  steady  flow 
process  will  be  a  function  of  only  the  upstream  reflected 
stagnation  conditions  as  long  as  the  dump  chamber  has 
sufficient  volume.   Ideally,  uniform  reservoir  stagnation 
conditions  will  end  when  the  reflected  region  conditions  are 
changed  by  a  wave  emanating  from  the  interaction  of  the 
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contact  surface  and  the  head  or  tail  of  the  reflected  expan- 
sion wave  (Figure  2).   The  useful  test  time  in  the  nozzle 
ends  when  the  reflected  pressure  drops  below  the  level 
required  to  maintain  the  static  pressure  at  the  nozzle  exit 
above  the  back  pressure  in  the  dump  chamber.   Thus  the 
useful  test  time  will  depend  on  the  nozzle  throat  size, 
which  governs  the  flow  rate  out  of  the  tube,  and  the  size  of 
the  dump  chamber,  which  determines  the  rate  at  which  the 
back  pressure  builds  up. 

D.   TAILORED  OPERATION 
1.   Tailoring; 

The  fundamental  limitation  of  a  shock  tube  is  the 
available  test  time.   Considerable  effort  has  gone  into 
trying  to  extend  shock  tube  test  times  [Ref.  10].   Test 
times  can  be  lengthened  by  as  much  as  a  factor  of  ten  by 
"tailoring,"  which  means  adjusting  conditions  so  that  the 
ratio  of  the  speed  of  sound  and  specific  heat  behind  the 
reflected  shock  matches  those  in  the  region  behind  the 
contact  surface.   When  this  tailored  condition  is  present, 
the  reflected  shock  passes  through  the  contact  surface 
bringing  it  to  rest.   A  Mach  wave  is  reflected  forwards  off 
the  contact  surface  and  generates  no  additional  disturbances 
as  it  passes  through  region  5.   In  the  untailored  condition, 
the  wave  reflected  forward  from  the  contact  surface  is  not  a 
Mach  wave  but  a  rarefaction  wave  when  ar   <    a3,  and  a  shock 
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wave  when  ax  >  a3.   The  reflected  waves  cause  conditions  in 
region  5  to  change.   Experimental  evidence  of  the  effect  of 
the  wave  interactions  obtained  by  Pennelegion  [Ref.  10]  are 
shown  in  Figure  10.   Pennelegion ' s  research  fully  explored 
tailoring  analytically  and  provided  supporting  experimental 
evidence.   His  findings  and  experimental  technigues  were 
followed  in  the  experiments  undertaken  in  the  present  study. 
The  tailored  interface  technigue  reguires  that  the  following 
eguation  be  satisfied  [Ref.lOrp.  4]; 


!i  =  !i 

a3   H 


y  +1   P     ■ 
2  Y]_  P2 


V1  P5 
1+  (2V)  (p^-1) 


L/2 

(3) 


A  unigue  tailored  condition  (shock  Mach  number) ,  exists  for 
every  combination  of  driver  and  driven  gas.   The  duration  of 
constant  pressure  in  region  5  lasts  until  the  rarefaction 
wave  reflects  from  the  driver  endplate  and  returns  to  meet 
the  contact  surface.   For  low  incident  shock  strengths  it  is 
the  head  of  the  rarefaction  wave,  and  for  high  incident 
shock  strengths  it  is  the  tail  of  the  rarefaction  wave  which 
terminates  the  constant  pressure. 

The  tailored  interface  technigue  for  producing 
extended  periods  of  constant  reflected  pressure  and 
temperature  was  demonstrated  experimentally  for  Mach  numbers 
below  3.4.   When  the  technigue  was  applied  to  shocks  of  Mach 
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number  greater  than  3.8,  there  was  a  significant  drop  in 
test  time  due  to  boundary  layer  interaction.   In  the  under- 
tailored  condition  the  reflected  conditions  showed  an 
initial  fall  in  temperature  and  pressure,  followed  by  a  slow 
rise  in  pressure  and  continued  drop  in  temperature.   The 
overtailored  condition  revealed  a  pressure  rise  from  the 
reflection  from  the  contact  surface  followed  by  a  rapid  drop 
in  temperature.   In  all  cases  with  strong  incident  shock 
waves,  the  real  gas  effects  caused  significant  variation  in 
the  reflected  region  conditions.   [Ref.  10] 
2 .   Mixing  of  Driver  Gases  (He  and  N2) 

The  application  of  tailoring  to  increase  test  time 
in  a  hypersonic  shock  tunnel  allows  more  accurate 
measurements  of  force  and  pressure  to  be  made  on  test 
models.   As  desired  test  section  Mach  number  increases  there 
is  a  significant  decrease  in  working  section  static 
pressure.   It  is  therefore  desirable  to  have  reflected 
region  pressures  as  steady  and  as  high  as  possible.  Pressure 
recovery,  defined  as  the  ratio  of  the  initial  step  in  the 
reflected  pressure  (P5)  to  driven  pressure  (PJ  ,  may  be 
increased  by  increasing  the  speed  of  sound  of  the  driver  if 
some  test  time  is  sacrificed.   The  combination  of  a  driver 
of  Helium  (He)  and  Nitrogen  (N2)  can  provide  an  improvement 
in  test  time  and  pressure  recovery  over  a  driver  of  N2 
alone.   Experimental  evidence  of  this  is  shown  in  Figures  11 
and  12  from  Pennelegion  [Ref.  10].   From  the  same  study, 
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Figure  13  and  Table  2  show  tailoring  using  various  mixtures 
of  He  and  N2  to  generate  shocks  in  the  range  1  <  Ms  <  3.4. 
The  tailoring  shock  Mach  number  may  be  selected  to  provide 
the  reguired  stagnation  temperature  to  prevent  liguefaction 
and  maximize  test  time  for  hypersonic  nozzles  with  various 
Mach  numbers.   Pennelgion  [Ref.  10],  provided  the  method  for 
calculating  tailored  conditions  for  He,  N2  mixtures. 
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TABLE  2 

EXPERIMENTALLY  TAILORED  MACH  NUMBERS  FOR  DRIVER 
GAS  MIXTURES  [REF.  10] 


He 

Time 

Steady 

T5 

Fract 

ion 

MTE 

P4/P1 

Pressure  (msec) 

(°K) 

0.7 

2.0 

24 

6.5 

714 

0.8 

2.7 

60 

6.0 

1096 

0.9 

2.8 

70 

4.7 

1159 

1.0 

3.4 

124 

4.3 

1553 
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III.   SHOCK  TUBE 

In  this  section  the  existing  shock  tube  and  instrumenta- 
tion, and  the  modifications  made  during  the  course  of  the 
present  experiments,  are  described. 

A.  STRUCTURE 

The  three  inch  diameter  shock  tube  at  the  Naval 
Postgraduate  School  Aeropropulsion  Laboratories  is  shown  in 
Figure  14.   A  detailed  description  is  given  in  Penaranda 
[Ref.  4].   The  tube,  constructed  of  highly  polished  303 
stainless  steel  consists  of  a  driver  section  7.75  feet  long, 
shown  in  detail  in  Figure  15,  and  a  driven  section  18.9  feet 
long. 

B.  DIAPHRAGM  CHAMBER 

The  tube  incorporates  a  double  diaphragm  arrangement  to 
maintain  precise  and  repeatable  diaphragm  bursting  pressure 
ratios  (P^/Px) .   The  tube  currently  is  fitted  with  diaphragms 
designed  to  burst  using  a  constant  driver  pressure  of  600 
psi.   Scored  brass  diaphragms  designed  to  rupture  at  400  psi 
are  inserted.  Variation  in  shock  speed  is  obtained  through 
adjustment  of  the  downstream  pressure  (Px)  .   To  obtain  the 
desired  pressure  ratio,  the  2.5  inch  long  diaphragm  chamber 
(Figures  15  and  16)  and  the  driver  chamber  are  simul- 
taneously filled  to  one  half  of  the  desired  P4  pressure  (300 

26 


Driver    Section- 


7.75' 


I    L 


-rront    cnfl 


-rront    End    Plote 


-Driven     Section 


Tube 


L 


Valve   ooay  Section 
Diapnrogm     Chamber 


18.9' 


61.375"    attd — 56.5 


Tube 


Tuce 


Instrumental  ion 


Stations 
A  B 


Pressure   Transducer 


'$ 


m 


'      7  I    Glow     oi 

P,ie<    Roil 

Tliin  Plallnum  Film 


A^P 


1 

A 

.     J 

1 

B 

Figure   14. 


3"    Diameter  Shock  Tube  with 
Instrumentation   Layout 


27 


Vacuum     Shul-off 
Valve 


High     Pressure     Gauge 


II   ~r   II 


High    Pressure       T       Connection 
Pressure     Control    Valve 


^— Quick    Acting    Valves 
-\\    (Connection    for 
two    bottles) 


Vacuum     Gauge 


Vacuum    Pump 


-^ 


Bleed  Valve 


Pressure *~ 

Regulator    Vnlve 


& 


Q 


(Q)    Front       or      Axial      View 


-*-High 
Pressure 
Gas 
Bottles 


Pressure     Gauge 

High   Pressure 
Control     Valve 

To  Diaphragm 
Chamber 


Flexible     Hose 


Driver     Section 


(b)    Side      View 


Bleed 
Valve 


Diaphragm 
Chamber 
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Figure  16.   Dual  Diaphragm  Chamber  (Unbolted) 

psi) .   A  diverting  valve  is  closed  to  isolate  the  diaphragm 
chamber  at  300  psi.   The  driver  section  pressure  is  then 
raised  to  the  desired  P^  pressure  (600  psi)  .   To  fire  the 
tube,  a  diverting  valve  is  opened  and  some  of  the  driver  gas 
at  600  psi  is  dumped  into  the  diaphragm  chamber  causing  the 
first  diaphragm  to  break.   The  second  diaphragm  breaks  im- 
mediately, as  the  pressure  on  the  driven  face  of  the 
diaphragm  rapidly  falls. 
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The  dual  diaphragm  arrangement  minimizes  bursting 
pressure  fluctuations  due  to  inconsistencies  in  the  material 
of  the  diaphragms.   By  applying  a  pressure  of  2  00  psi  above 
the  design  bursting  pressure  of  the  first  diaphragm,  rupture 
of  the  two  diaphragms  is  ensured  with  a  precise  value  of  PA. 
In  shock  tubes  using  a  single  diaphragm  arrangement,  precise 
bursting  pressures  cannot  be  repeated  unless  the  complexity 
of  a  diaphragm  puncturing  device  is  used.   Without  a 
puncture  device,  the  bursting  pressure  fluctuations  remain  a 
function  of  the  uniformity  of  the  diaphragm  material. 

C.   NOZZLE  ASSEMBLY 

An  existing  two  dimensional  stainless  steel  nozzle 
designed  for  an  exit  Mach  number  of  5.72  (A/A*  =  4  3.7)  and  a 
rectangular  dump  chamber  to  which  it  bolted  were  available 
for  the  present  study.   The  assembly  is  shown  in  Figure  17. 
For  the  present  experiments  two  1"  diameter  window  ports 
were  installed  to  allow  optical  measurement  techniques  to  be 
exercised,  and  a  model  support  was  constructed  to  position  a 
3/8"  axisymmetric  model  on  the  nozzle  centerline  in  line 
with  the  window,  where  the  Mach  number  was  about  4.3.   The 
final  assembly  is  shown  in  Figure  18.   The  dimensions  of  the 
nozzle  are  given  in  Figure  19  and  the  design  of  the  window 
installation  is  shown  in  Figure  20.   Views  of  the  nozzle  and 
model  support  are  shown  in  Figure  21.   Since  the  nozzle  and 
dump  tank  existed  (from  a  previous  gas  dynamic  laser 
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Figure    18.      Assembled   Dump   Chamber,    Model    Support 
and   Nozzle 
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Figure  19.   Dimensional  Drawing  of  the  Nozzle 
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(a)   Model,  Model  Support  and  Nozzle 


(b)   Nozzle  Exit 


Figure  21.   Views  of  the  Nozzle  and  Model  Support 
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investigation) ,  location  and  sizing  of  the  window  as  well  as 
the  maximum  flow  duration  were  limited  by  the  existing 
hardware. 

D.   INSTRUMENTATION 
1.   Shock  Speed 

Shock  speed  was  determined  by  measuring  the  time  for 
the  shock  to  travel  between  thin  film  gauges  located  at 
positions  A  and  B  and  from  B  to  the  pressure  transducer   at 
C  (Figure  14) .   The  locally  manufactured  thin  film  gauges 
shown  in  Figure  22,  consist  of  a  0.25"  diameter  pyrex  rod 
coated  with  a  1/32"  wide  thin  film  of  platinum  giving  a  cold 
resistance  of  100  ohms.   [Ref.  4:p.  23]   As  the  region  of 
hot  gas  immediately  behind  the  shock  front  passes  over  the 
platinum,  its  temperature  and  resistance  change  and  a 
voltage  pulse  is  generated.   In  the  present  work  the  signals 
from  the  thin  films  were  amplified  and  recorded  for 
analysis.   The  thin  films  at  location  A  also  served  as  a 
trigger  for  the  digital  counter  and  shadowgraph  spark 
source. 

Shock  speed  was  determined  in  different  ways.   The 
pre-existing  techniques  involved  using  oscilloscope  records 
of  thin  film  and  pressure  transducer  outputs,  and  also  the 
output  of  a  digital  counter  triggered  on  and  off  by  the  thin 
films  at  A  to  B  respectively.   In  the  present  work,  in  runs 
1  through  49  shock  speeds  were  determined  from  a  Polaroid 


35 


Shock  Tube  Instrumentation 
Station  Body 
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Instrumentation 
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Figure  22. 


Thin  Film  Instrumentation  Assembly 
[Ref.  4] 
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photograph  of  the  oscilloscope,  such  as  is  shown  in  Figure 
23,  and  from  the  digital  counter.   Both  techniques  were 
identical  to  those  used  in  the  late  1960's. 

In  runs  50-8  5,  the  outputs  of  the  three  traces  were 
recorded  digitally  on  three  of  the  available  four  channels 
of  a  Scientific  Alanta  SD380  Spectrum  Analyzer.   Data  were 
recorded  at  .019  and  .039  msec  intervals  depending  upon  the 
time  scale  selected.   An  example  is  shown  in  Figure  24.   The 
shock  speeds  were  determined  during  postprocessing  of  the 
data . 

2 .  Driver  Tube  Pressure 

A  dial  gauge  with  10  psi  pressure  increments  was 
used  to  determine  when  the  constant  600  psi  driver  pressure 
was  reached.   A  dead  weight  tester  was  used  to  calibrate  the 
gauge . 

3 .  Driven  Tube  Pressure 

Two  calibrated  Wallace  and  Tiernan  Bourdon  tube 
gauges  were  used  to  determine  the  driven  tube  pressures.   In 
the  range  from  0  to  50mm  Hg,  a  gauge  with  increments  of  .2 
mm  Hg  was  used.   For  higher  driven  pressures,  a  second  gauge 
with  a  range  from  0  to  35  psi  in  increments  of  .05  psi  was 
used . 

4 .  Driven  Tube  Temperature 

A  Type  K  (chromel-alumel )  thermocouple,  connected  to 
a  Daigger  Scientific  Digital  Panel  Meter  using  matched 
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Figure  23. 


Typical  Polaroid  Photograph  of  the 
Thin  Film  Response  from  Location  B 
and  the  Output  of  the  Original 
Pressure  Transducer 
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extension  wires,  was  mounted  in  the  driven  end  of  the  tube. 
The  digital  readout  was  recorded  prior  to  every  run. 

5 .  Reflected  Pressure  and  Test  Duration 

A  piezoelectric  pressure  transducer,  Endevco  model 
2501,  a  charge  amplifier,  Endevco  model  2700,  and  a 
Tektronics  551  oscilloscope  were  initially  used  to  record 
the  transient  region  5  pressure.   The  lack  of  resolution  and 
strange  gualitative  behavior  observed  in  the  results  (Figure 
23),  led  to  the  mounting  of  an  improved  transducer. 

A  PCB  Piezontonics,  Inc.  Model  124A  pressure  trans- 
ducer connected  to  a  model  504A  dual  gain  charge  amplifier 
were  used  in  runs  50  through  75.   The  system  was  calibrated 
using  a  standard  static  calibration  feature  described  in  the 
manufacturer's  manuals  [Refs.  11,12].   The  transducer 
provided  a  response  of  240  psi/volts.   Experimental  traces 
obtained  are  given  in  Appendix  D. 

6 .  Oscilloscopes 

Two  Tektronics  oscilloscopes,  models  549  and  551, 
fitted  with  Polaroid  cameras  were  used  for  runs  1  through 
49.   The  capability  of  the  scopes  are  fully  described  in  the 
associated  manuals  [Refs.  13,14].   Time  measurements  to 
determine  shock  speed,  and  voltage  levels  to  determine 
pressure,  were  determined  by  overlaying  a  scale  on  the 
Polaroid  photographs. 
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7 .   Signal  Analyzer 

A  Scientific  Alanta  SD380  four  channel  signal 
analyzer  (Figure  25)  equipped  with  a  hard  disk  and  plotter 
was  used  to  record  and  process  data  from  runs  50  through  85. 
The  independent  channel  trigger  delays  (±) ,  adjustable 
scales  and  post-processing  capabilities  of  the  digitized 
data  (in  contrast  to  photographing  oscilloscopes),  is 
evidence  of  the  advances  in  measurement  capabilities  over 
the  past  20  years.   A  full  description  of  capabilities  of 
the  analyzer  is  given  in  the  instrument  manual.   [Ref.  15] 

E.   OPTICAL  MEASUREMENT  SYSTEMS 

A  simple  shadowgraph  arrangement  was  generated  from  the 
components  of  a  6"  schlieren  system.   The  optical  layout  is 
shown  in  Figure  26  [Ref.  16].   The  spark  source  was  used 
with  Polaroid  film  in  a  camera  back  attached  to  the  nozzle 
window.   The  continuous  source  was  used  with  a  Hycam  movie 
camera. 

1 .   Light  Sources 

a.   High  Pressure  Mercury  Arc  Lamp 

A  high  pressure  mercury  arc  lamp  (Per  Type  112- 
5145)  constant  light  source  was  used  with  a  six  inch 
achromatic  convex  lens  with  a  focal  length  of  60  inches.  The 
arrangement,  mounted  on  an  adjustable  bench  (Figure  27) , 
provided  a  parallel  light  beam  through  the  nozzle  windows. 
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Figure  25.  Scientific  Alanta  SD380  Four 
Channel  Signal  Analyzer  with 
Hard  Disk  Installed 
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Figure  26.   Optical  Layout  Sketch 
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Figure  27.   Constant  Light  Source  and  Lens  Mounted 
on  an  Adjustable  Optical  Bench 


Figure  28. 


Spark  Source  with  High  Voltage 
Power  Supply 
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The  beam  was  adjusted  by  the  procedures  described  in  the 
associated  manual  [Ref.  16]. 

b.   NASA  229-D  Spark  Light  Source 

A  NASA  229-D  spark  light  source  (Figure  28) 
provided  a  nanosecond  spark  for  shadowgraph  pictures  taken 
with  Polaroid  film.   The  spark  was  triggered  from  a  thin 
film  connected  through  the  delay  and  voltage  circuits  shown 
in  Figure  29.   The  timing  for  the  spark  source  was 
determined  experimentally. 
2 .   Camera  Systems 

a.  Hycam  Movie  Camera 

A  Hycam  movie  camera  [Ref.  17]  with  a  framing 
rate  on  the  order  of  6000  pictures  per  second  at  runout  of  a 
400  ft  roll  of  film,  was  used.   The  camera  was  equipped  with 
a  Nikon  55mm  lens  and  ASA  250  black  and  white  film.   A  time 
reference  was  provided  on  the  edge  of  the  film  by  a  strobe 
light  pulsing  at  1000  pulses  per  second.   The  camera  was 
triggered  manually.   The  camera  was  mounted  on  a  second 
adjustable  bench  as  shown  in  Figure  30. 

b.  Spark  Shadowgraphs 

Spark  shadowgraphs  were  taken  using  a  Polaroid 
Land  Film  Holder  with  Type  52  ASA  400  Black  and  White  film, 
held  directly  behind  the  optical  windows.   To  eliminate 
stray  light,  the  mounting  arrangement  was  contained  in  a 
cardboard  box,  and  tests  were  conducted  after  dark. 
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Figure    29. 


Spark  Timing  Delay  Circuit  and  High 
Voltage  Pulse  Generator 


Figure  30.   Hycam  Movie  Camera  with  Separate 
Pulse  Generator 
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IV.   EXPERIMENTAL  PROCEDURE 

A.  PREPARATION 

Prior  to  testing,  with  the  driven  endplate  removed  and 
the  diaphragm  chamber  opened,  a  cloth  lightly  dampened  with 
Acetone  was  pulled  through  the  driver  and  driven  sections. 
O-ring  grooves  of  the  opened  sections  where  thoroughly 
cleaned  and  a  thin  coating  of  vacuum  grease  was  applied  to 
the  rubber  O-rings.   Diaphragms  were  wiped  clean  of  visible 
dirt  and  the  system  reassembled.   The  cleaning  procedure  was 
repeated  after  approximately  15  runs.   Using  a  roughing 
pump,  a  clean-tube  vacuum  pressure  of  3  mm  Hg  was  routinely 
achieved. 

B.  LOADING  THE  DRIVEN  GAS 

Prior  to  loading  to  the  desired  test  pressure  (Px)  ,  the 
driven  section  was  purged.   To  purge,  the  tube  was  evacuated 
to  3  mm  Hg  and  then  isolated  from  the  roughing  pump. 
Nitrogen  was  introduced  to  a  pressure  1  to  4  psi.   The  tube 
was  then  evacuated  a  second  time  and  repressurized  to  the 
desired  test  condition. 

C.  LOADING  THE  DRIVER  GAS  AND  FIRING 

The  driver  assembly  is  designed  to  operate  at  600  psi 
with  single  or  dual  gas  combinations.   To  purge  the  driver 
section  and  diaphragm  chamber,  the  system  was  pressurized  to 
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approximately  30  psi  and  then  vented  to  atmospheric  pres- 
sure.  The  system  was  then  pressurized  to  600  psi  with  the 
diverter  valve  to  the  diaphragm  chamber  being  closed  at  300 
psi.   When  the  70%  He  and  30%  N2  driver  mixture  was  used, 
the  system  was  pressurized  to  420  psi  with  He  and  then  to 
600  psi  with  N2. 

The  tube  was  fired  by  opening  the  diverter  valve  to  the 
diaphragm  chamber.   Some  seepage  from  the  driver  section  was 
evident.   To  minimize  this  effect  the  tube  was  fired 
immediately  after  reaching  the  600  psi  level.   Replacement 
of  O-rings  would  minimize  this  problem.   Seepage  into  the 
driven  section  was  not  evident. 
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V.   EXPERIMENTAL  RESULTS  AND  DISCUSSION 

A  series  of  experiments  was  conducted  to  gain  experience 
with  the  shock  tube.   The  study  consisted  of  approximately 
85  shots  involving  the  use  of  He  and  a  70%  He/30%  N2  mixture 
as  driver  gases,  driven  into  a  N2  test  gas.   Studies  with 
the  shock  tube  assembled  with  the  pressure  transducer 
mounted  in  the  end  of  the  driven  tube,  were  concentrated  on 
investigating  the  concept  of  tailoring  to  maximize  the 
magnitude  and  duration  of  steady  pressure  in  the  reflected 
region.   While  more  than  75  shots  were  made  in  this 
configuration,  the  replacement  of  the  pressure  transducer 
and  postprocessing  of  data  taken  with  the  SD380  greatly 
enhanced  the  results.  For  this  reason,  only  data  taken  after 
run  50  using  the  SD380  are  presented.  Studies  with  the 
nozzle  attached  at  the  end  of  the  driven  tube  were 
exploratory  attempts  to  obtain  flow  visualization  in  a  shock 
tunnel  configuration.  The  results  of  the  two  series  of 
experiments  are  discussed  in  the  following  sections. 

A.   TAILORING 

Using  the  nomenclature  of  Glass  [Ref.  18],  a  tailored 
condition  is  defined  as  the  condition  when  a  Mach  wave  is 
reflected  off  the  contact  surface  and  causes  no  change  in 
the  reflected  pressure  (P5)  .  This  can  be  identified 
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experimentally  (Figure  10)  as  the  condition  when  the  ratio 
of  the  initial  reflected  pressure  (P5)  and  the  pressure 
after  the  wave  interaction  at  the  contact  surface  reflects 
off  the  endplate  (P8)  being  egual  to  one.   Figure  31  shows 
examples  of  results  in  the  present  study  for  the  end  wall 
pressure  untailored  and  tailored  conditions  using  He  with  an 
incident  shock  strength,  Ms  =  3.4.   Similar  traces  for  near 
to  a  tailored  mixture  (.7  He  and  .3  N2)  with  Ms  =  2.0  are 
shown  in  Figure  32.   These  results  are  in  complete  agreement 
with  Pennelegion  [Ref.  10]. 

B.   TIME  OF  STEADY  PRESSURE 

It  was  experimentally  determined  in  the  present  study 
that  a  tailored  mixture  as  a  driver  gas  produced  a  longer 
duration  of  steady  pressure  than  that  of  a  pure  He  driver. 
The  time  of  steady  pressure  for  the  existing  NPS  shock  tube 
using  He  as  a  driver  was  approximately  3.6  milliseconds 
while  the  mixture  extended  the  test  time  duration  to  5 
milliseconds.   A  comparisons  of  the  test  times  obtained  in 
the  present  study  with  those  of  Pennelegion  [Ref.  10]  is 
shown  in  Figure  33.   The  trends  are  seen  to  be  similar. 
Since  the  work  of  Pennelegion  [Ref.  10]  was  performed  at  a 
constant  driver  pressure  of  1000  psi  in  a  N.P.L.  2"  shock 
tube  with  a  10  ft  driver  and  a  12  foot  driven  section  some 
differences  in  the  magnitudes  of  the  test  times  are  to  be 
expected. 


50 


V  ! 


t 


ur 


j^i 


•jjw^^'^Sw. 


REFLECTED  PRESSURE  RESPONSE 
HE  DRIVER/NITROGEN  DRIVEN 


-2 


RUN  58   (MACH  =  3.8) 


(a)   Untailored 


J 


AT  .019m  SEC    .010000 


r 


fVfc, 


^V^\^^VVvvv^M^M^WvVV^ 


n 

REFLECTED  PRESSURE  RESPONSE 
HE  DRIVER/NITROGEN  DRIVEN 

RUN  60   (MACH  =3.8) 

i 

2 

t 

i 

( 

AT  .019m  SEC 

.010000 

(b)   Tailored 


Figure  31. 


Experimentally  Determined  Endwall  Pressure 
Traces  Using  Helium  as  a  Driver  Gas 
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Figure  32. 


Experimentally  Determined  Endwall  Pressure 
Traces  Using  a  Mixture  of  70%  Helium  and 
30%  Nitrogen  as  a  Driver  Gas 
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For  low  strength  shocks  the  head  of  the  expansion  wave 
causes  an  early  degradation  of  the  reflected  pressure  in  the 
NPS  tube.   The  longer  driver  section  in  the  N.P.L.  tube  is 
responsible  for  the  difference. 

C.  PRESSURE  RECOVERY 

Figure  34  shows  the  fraction  of  initial  driver  pressure 
available  in  the  reflected  region  as  a  function  of  shock 
Mach  number.   The  data  show  that  as  shock  strength  increases 
pressure  recovery  decreases  for  both  He  and  mixed  driver 
gases.   The  data  are  in  good  agreement  with  the  results  of 
Pennelegion  [Ref .  10] .   The  measured  values  of  the  reflected 
shock  pressure,  which  are  the  pressures  available  to  drive  a 
hypersonic  tunnel,  are  given  in  Tables  3  and  4. 

D.  REFLECTED  TEMPERATURE 

The  reflected  temperature  was  not  determined 
experimentally.   Instead,  shown  in  Table  2  are  values  of 
stagnation  temperatures  calculated  for  shock  strengths  in 
the  range  from  2.0  to  3.4.   These  values  were  taken  from  the 
tabulated  solutions  (based  on  equilibrium  chemistry)  of 
Bernstein  [Ref.  19:pp.  32-43]. 

E.  SHOCK  ATTENUATION 

A  representative  illustration  of  measured  shock 
attenuation  for  both  a  pure  gas  driver  and  the  driver  gas 
mixture  is  shown  in  Figure  35.   Velocity  can  be  inferred 
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Figure  34.   Comparison  of  Reflected  Pressure  Measurements 
in  the  NPS  3  Inch  Shock  Tube  with  Experimen- 
tal Results  of  Pennelegion  [Ref.  10] 
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TABLE  3 

EXPERIMENTALLY  DETERMINED  REFLECTED  PRESSURE  MAGNITUDE 
AND  DURATION  AND  CALCULATED  TEMPERATURE  (HE  DRIVER/N2  DRIVEN) 


un 

P5 

Time  of  Steady 

# 

Ms 

P,/Pi 

(psi) 

P8/P5 

P5  (msec) 

50 

3.08 

994 

150 

1.79 

1.3 

51 

2.92 

353 

235 

1.23 

1.2 

52 

4.39 

414 

207 

1.38 

1.9 

53 

3.32 

93 

312 

1.42 

1.5 

57 

4.52 

500 

183 

1.25 

0.3 

58 

3.79 

207 

304 

1.23 

0.3 

59 

3.66 

148 

386 

1.13 

1.6 

60 

3.36 

169 

264 

1.13 

3.6 

61 

3.72 

167 

317 

1.23 

0.5 

64 

3.42 

164 

— 

65 

3.65 

158 

— 

66 

184 

67 

3.77 

173 

68 

3.81 

400 

74 

4.50 

667 

75 


4.10* 


1240 


*  Estimated  Value 


TABLE  4 

EXPERIMENTALLY  DETERMINED  PRESSURE  MAGNITUDE  AND  DURATION 
AND  CALCULATED  TEMPERATURE  (70%  HE,  30%  N2  DRIVER/N2  DRIVEN) 


Run 

P5 

Time 

of  Steady 

# 

Ms 

P4/P1 

(psi) 

P8/P5 

P, 

5  (msec) 

62 

2.40 

50 

317 

1.10 

2.8 

63 

2.66 

51 

310 

1.23 

1.5 

69 

2.47 

50 

310 

1.09 

2.5 

70 

2.25 

50 

305 

1.11 

3.3 

71 

2.03 

24 

427 

1.02 

4.9 

72 

2.03 

22 

422 

1.08 

5 

73 

1.99 

21 

410 

1.13 

2 

77 

78 

155 

79 

3.0* 

171 

80 


3.4* 


638 


*  Estimated  Value 
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from  the  slope  of  the  lines  between  the  thin  films  A  and  B 
and  between  B  and  the  pressure  transducer  at  C.   An 
increasing  slope  indicates  a  slower  speed.   The  general 
trend  is  that  for  a  given  gas  mixture,  the  stronger  the 
shock  strength  the  smaller  the  attenuation  between  B  and  C. 
Reflected  velocity  is  of  the  order  of  one  half  the  incident 
velocity.   Tabulated  results  for  all  runs  are  found  in 
Appendix  D. 

F.   OPTICAL  MEASUREMENTS 

Before  attempting  to  obtain  photographs  during  the 
firing  of  the  shock  tube,  a  small  freejet  from  a  1/4" 
flattened  copper  tube,  powered  by  shop  air,  was  installed  in 
the  model  location,  the  spark  light  was  adjusted  using  a 
neutral  filter  and  a  successful  shadowgraph  was  taken. 

1 .   Shadowgraph  Using  Hycam  Camera 

The  Hycam  Movie  camera  was  triggered  manually  either 
one  or  two  seconds  before  the  diaphragm  was  ruptured.   The 
400  ft  rolls  of  ASA  250  speed  film  provided  approximately 
one  second  for  the  camera  to  accelerate,  and  a  three  second 
window  to  capture  the  event.   It  was  found  that  all  the 
available  light  was  needed  to  obtain  a  correct  exposure  and 
the  camera  was  focused  on  the  outside  surface  of  the  window. 
The  shadowgraph  was  made  visible  by  the  faint  film  which  by 
then  had  deposited  on  the  inside  of  the  window.   The  result 
was  successful. 
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With  the  blunt  body  shown  in  Figure  21(a)  installed, 
the  nozzle  starting  Mach  wave  (Figure  36(a)  and  (b)  )  ,  and 
the  establishment  of  steady  flow  with  an  oblique  shock 
visible  (Figure  36(c)),  were  successfully  recorded.   Steady 
flow  (started)  conditions  lasted  a  brief  time  until  dump 
chamber  back  pressure  unstarted  the  nozzle.   Incident  shock 
strength  was  estimated  to  be  Ms  =  3 . 8  providing  a  peak 
pressure  in  the  reflected  region  of  approximately  460  psi. 
The  actual  shock  speed  was  not  measured  during  the  test 
since  the  starting  of  the  camera  caused  premature  triggering 
of  the  spectrum  analyzer.   However  the  initial  conditions, 
and  therefore  the  reflected  pressure,  were  the  same  as  in 
run  58.   A  complete  record  of  this  run  given  in  Appendix  D. 
2 .   Spark  Shadowgraph 

Before  firing  the  shock  tube,  the  correct  exposure 
for  the  expected  shadowgraph  was  verified  by  opening  the 
film  and  firing  the  spark  source.   The  result  is  shown  in 
Figure  37.   It  was  subsequently  found  that  when  the  shock 
tube  was  fired,  the  thin  film  signal  failed  to  trigger  the 
spark,  and  thus  any  exposure  on  the  film  was  the  result  of 
light  produced  from  within  the  tube.   The  results  of 
exposing  the  film  without  discharging  the  spark  as  the  shock 
Mach  number  was  progressively  increased,  are  shown  in 
sections  (a)  to  (e)  of  Figure  38.   A  faint  glow  appeared 
first  at  a  ME  =  3.6.   As  Ms  was  increased  for  both  a  pure  He 
driver  and  mixed  driver  gases,  the  luminosity  increased. 
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Nozzle  Starting  Shock  Wave 


(a)      Nozzle  Starting  Shock  Wave 


(b)     Nozzle  Starting  Shock  Wave 

/<-  Oblique  Shock  Wave 


(c)      Oblique  Shock 


Figure    36.      Results   Obtained  Using  Hycam   Camera 
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Figure  37.   Spark  Shadowgraph  with  No  Flow  Present 
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The  source  of  the  light  requires  further  investigation,  but 
it  is  clear  that  some  contamination  must  be  present  since 
the  driver  and  driven  gases,  if  pure,  would  not  be  excited  at 
the  temperature  produced  in  the  firing  process. 
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VI.   SIZING  CONSIDERATIONS  FOR  A  HYPERSONIC  TUNNEL 

The  experimental  program  showed  that  the  basic  shock 
tube  gave  reliable  and  repeatable  results  and,  following  the 
upgrading  recommendations  detailed  in  Appendix  A,  it  could 
well  serve  as  a  driver  for  a  hypersonic  tunnel.   Having 
experimentally  defined  the  stagnation  temperature  and 
pressure  available  to  drive  such  a  facility,  it  only 
remained  to  arrive  at  a  nozzle  and  a  dump  chamber 
configuration  suitably  matched  to  the  tube.   Simplified 
models  were  therefore  developed  to,  first,  calculate  the 
effect  of  nozzle  throat  diameter  on  the  available  steady 
flow  time  and,  second,  to  determine  a  suitable  dump  chamber 
size.   The  two  calculations  are  detailed  in  Appendix  B  and 
Appendix  C  respectively. 


A.   NOZZLE  THROAT  SIZE 

The  nozzle  throat  diameter  must  be  sized  such  that  the 
time  for  which  conditions  in  the  reflected  region  are  steady 
must  be  greater  than  the  time  required  for  the  compressed 
mass  of  gas  to  exhaust  through  the  nozzle.   The  model 
described  in  Appendix  B  estimates  the  time  in  relation  to 
the  falling  pressure  in  the  reflected  region  to  be  given  by 


t  = 


C1A*P5 


m 


Y5"l 


- 1 


2y« 


LP5 


(B.10) 
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If  it  is  assumed  the  supply  pressure  can  drop  to  half  its 
initial  value,  the  time  required  can  be  calculated  using 
Equation  (B.10)  for  different  driver  gases  and  shock 
strengths.   Through  the  ratio  m0/m0,  the  calculated  time  is 
implicitly  a  function  of  the  throat  diameter.   Results  of 
such  calculations  are  shown  in  Figure  39(a)  and  (b) .   The 
four  curves  shown  illustrate  the  dependence  of  useful  time 
on  the  shock  strength  and  throat  diameter.   Clearly,  the 
shock  strength  is  less  important  than  the  throat  area. 
Based  on  these  results,  a  throat  diameter  6.5mm  (.025 
inches)  was  selected.   This  choice  made  best  use  of  the 
steady  flow  time  which  was  experimentally  determined  to  be 
from  3.5  to  5  milliseconds  using  tailored  conditions.   For 
the  6.5mm  throat,  Equation  (B.10)  was  used  to  calculate  the 
time  for  the  reflected  shock  region  pressure  to  decay  to 
various  fractions  of  its  initial  magnitude.   The  results  are 
shown  in  Figure  40.   Clearly,  for  even  the  strongest  shock, 
the  time  for  the  pressure  to  decay  to  80%  of  its  initial 
value  is  longer  than  its  experimentally  determined  tailored 
test  time  of  5  msecs. 

B.   DUMP  CHAMBER  SIZE 

The  dump  chamber  requires  a  capacity  such  that  the  back 
pressure  remains  below  the  nozzle  exit  pressure  for  the 
duration  of  steady  conditions  in  the  driven  tube.   The 
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(a)   Helium  into  Nitrogen 


Figure  39.   Time  for  Reflected  Region  Pressure  to  Drop 
to  Half  its  Initial  Value 
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Figure  39  (CONTINUED) 
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Figure  40.   Calculated  Time  of  Decay  of  Reflected 

Region  Pressure  for  a  6.5mm  Nozzle  Throat 
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analysis  in  Appenix  C  gave  the  required  ratio  of  dump 
chamber  volume  to  reflected  shock  reservoir  volume  as 


V, 
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Examination  reveals  that  the  parameters  which  most  influence 
the  required  dump  chamber  volume  are  the  incident  shock 
strength  (which  determines  the  reflected  reservoir 
conditions) ,  the  desired  test  section  Mach  number,  the 
driver  gas  selected,  and  the  specified  reflected  region 
pressure  (PL/Pt0)  at  which  the  flow  breaks  down.   Since  the 
pressure  recovery  is  a  function  of  which  driver  gas  is 
selected,  calculations  were  made  for  both  nitrogen  and 
helium.  The  results  are  shown  in  Figure  41.   A  mixture  of 
the  driver  gases  would  be  expected  to  fall  between  these  two 
curves.   Consistent  with  the  finding  above,  the  criterion 
that  the  dump  chamber  pressure  becomes  equal  to  the  nozzle 
exit  pressure  when  20%  of  the  stagnation  pressure  was 
exhausted,  was  taken  to  determine  the  necessary  dump  chamber 
volume.   Shown  in  the  figure  are  results  obtained  assuming 
the  unknown  temperature  in  the  dump  chamber  to  have  its 
extreme  maximum  and  minimum  values;  namely,  initial 
reflected  shock  temperature  and  initial  dump  chamber 
temperature.   From  an  examination  of  the  results,  a 
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Figure  41.   Dump  Chamber  Volume  Requirements 
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conservative  selection  would  be  15  cubic  meters.   This  is  a 
vessel  approximately  6  feet  in  diameter  and  15  feet  long. 

C.   HIGH  ENTHALPY  MODIFICATION 

Higher  enthalpy  flows  can  be  generated  using  the  free- 
piston  arrangement  shown  in  Figure  1.   Stalker  [Refs.  20,21] 
has  carried  out  extensive  experimental  work  in  this  area.   A 
preliminary  analysis  shows  that  the  NPS  shock  tube  has  the 
same  internal  diameter  as  Stalker's  1st  generation  free- 
piston  tunnel,  but  has  a  thinner  wall  by  1/2".   Structural 
limitations  would  prevent  the  full  duplication  of  Stalker's 
tunnel  but,  with  the  published  information  on  the  design  of 
free  piston  drivers,  it  should  be  possible  to  extend  the 
present  tube  to  an  increased  enthalpy  level  at  which  real 
gas  effects  of  interest  can  be  studied. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 


The  present  study  showed  that: 

The  existing  shock-tube  is  straightforward  and 
relatively  inexpensive  to  operate,  and  gives 
repeatable  test  conditions; 

The  existing  shock  tube  can  be  used  to  drive  a 
hypersonic  shock  tunnel  to  a  Mach  number  of  12 ; 

Tailored  operation  can  provide  uniform  stagnation 
conditions  for  periods  of  approximately  3.4  msec 
using  Helium,  and  approximately  5  msec  using  a  .7 
He/.3N2  mixture,  as  the  driver  gas; 

It  is  feasible  to  obtain  both  high  speed  movies  and 
short  duration  spark  flow  visualization  records  of 
the  test  section  flow  field  using  currently  available 
equipment; 

Contamination  in  the  tube  can  lead  to  unwanted 
radiation  (which  must  clearly  be  eliminated  if  real 
gas  effects  are  the  purpose  of  the  experiment) . 

It  is  recommended  that: 

The  existing  tube  be  serviced  and  its  operating 
instrumentation  be  upgraded,  as  outlined  in  Appendix  A; 

A  Mach  10  nozzle,  15  cubic  meters  dump  chamber  and 
model  support  be  designed  and  procured  for  operation 
with  the  current  tube; 

Consideration  be  given  to  the  addition  of  a  free- 
piston  tunnel  to  provide  high  enthalpy  simulation. 

Operation  of  a  shock  tunnel  with  the  existing  tube  would 

provide  a  convenient  and  inexpensive  facility  for  use  in 

laboratory  instruction.   It  would  also  provide  a  tool  for 

thesis  studies  aimed  at  hypersonic  CFD  code  validation. 


74 


Extension  to  a  free-piston  driver  would  provide  an  advanced 
tool  for  hypersonic  research. 
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APPENDIX  A 
APPRAISAL  OF  THE  PRESENT  SHOCK  TUBE 

The  NPS  shock  tube  has  been  in  operation  since  the  late 
1960's.   This  appendix  provides  a  description  of  recommended 
improvements/enhancements  of  the  existing  system,  which 
result  from  the  experience  gained  in  the  reported  program  of 
tests. 

A.  O-RINGS 

The  planned  relocation  of  the  facility  to  the  new  Gas 
Dynamics  Laboratory  requires  complete  disassembly  of  the 
shock  tube.   Many  O-rings  are  old  and  worn  and  they  should 
be  replaced  at  that  time.   Table  5  details  the  O-ring 
requirements . 

B.  THIN  FILM  GAUGES 

The  unavailability  of  replacement  thin  film  gauges  used 
to  measure  shock  speed  is  unacceptable.   Only  the  two  in- 
stalled gauges  remain  operational.   The  options  are  either 
to  manufacture  thin  films  in-house  as  detailed  in  Penaranda 
[Ref.  4]  or  to  procure  them  commercially  through  the  supply 
system. 
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TABLE  5 
O-RING  REPLACEMENT  GUIDE 


Number 

Width  of 

Depth  of 

Inside 

Outside 

Location 

Required 

Groove 

Groove 

Diameter 

Diameter 

(inches) 

(inches) 

(inches) 

(inches) 

Driver  Tube 

1 

.100 

.157 

3.502+.005 

.106 

.16  3 

Space/ 

2 

.100 

.157 

4.000 

Diaphragm 

.106 

.16  3 

4.005 

2 

.100 

.157 

3.500 

.106 


.163 


3.505 


Valve  Body      1 

.101 

.157 

3.2485 

.105 

.163 

3.2495 

Tubing         3 
Sections 

.101 
.105 

.157 
.105 

3.2485 
3.2495 

Instrumentation  5 
Section 

.101 
.105 

.157 
.16  3 

3.2485 
3.2495 

Driver  End  Cap  1 

.101 

.157 

3.2470 

105 


16  3 


3.2460 
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C.  PRESSURE  TRANSDUCERS 

1.  Piezoelectric  Transducer 

The  two  Endevco  Model  2  501  transducers  with  Model 
2642  M35  charge  amplifiers  should  be  returned  to  the  manu- 
facturer for  calibration.   Biased  and  unbiased  responses 
should  be  checked.   Unknown,  but  apparently  unrealistic 
response  characteristics  made  these  transducers  unusable  for 
the  reported  experiments. 

2 .  Piezoresistive  Pressure  Transducer 

Procurement  of  a  piezoresistive  pressure  measurement 
system,  including  a  transducer  and  amplifier  such  as  the 
Kistler  Model  603B  Series  and  Model  5004  Charge  Amplifier, 
would  provided  continuous  measurements  of  shock  tube 
pressures. 

A  second  transducer,  to  be  used  as  a  pitot  probe  for 
flow  survey  work  in  shock  tunnel  applications,  is  also 
required.   Selection  would  be  based  on  the  selected 
nozzle  Mach  number  and  probe  geometry.   Cornell  Aeronautical 
Laboratory  publications  by  MacArther  should  be  consulted. 

D.  THIN  FILM  SIGNAL  AMPLIFIER 

The  current  thin  film  signal  amplifier  has  no  record  of 
calibration  for  the  past  ten  years.   Voltage  fluctuations 
are  apparent  and  dial  adjustment  response  is  very  poor. 
Prior  to  the  next  use  the  amplifier  components  should  be 
tested  and  replaced  as  required. 
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E.  DATA  ACQUISITION 

It  is  recommended  that  the  Scientific  Alanta  signal 
analyzer  or  its  equivalent  be  used  in  all  future  testing. 
Post  processing  capability,  and  the  availability  of  four 
channels  with  independent  trigger  adjustment  provide  a 
valuable  experimental  capability  not  available  with  the 
Tektronics  549  and  551  oscilloscopes. 

F .  DIAPHRAGMS 

The  existing  brass  diaphragms  allow  driver  operation  at 
a  constant  600  psi  level.   Redesign  to  a  higher  operational 
level  of  approximately  700  psi  would  increase  the  range  of 
reflected  pressures  available  to  drive  a  shock  tunnel,  with 
an  associated  cost  of  higher  gas  consumption.   While  not 
essential,  this  modification  would  reduce  problems  which 
would  be  expected  in  an  extension  to  a  shock  tunnel  with 
respect  to  the  measurement  of  small  aerodynamic  forces. 

Gaydon  [Ref.  2:p.  91]  and  Glass  [Ref.  18:  pp.  503-511] 
provide  valuable  design  insight  into  successful  diaphragm 
material  selection  and  preparation. 

G.  MISCELLANEOUS  HARDWARE 

Valve  7  and  valve  11  located  in  the  driven  end  of  the 
tube  [Ref.  4 : Figure  2]  exhibited  some  signs  of  leakage 
during  the  loading  process.   Repair  or  replacement  of  the 
valve  seats  is  recommended. 
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APPENDIX  B 
ESTIMATION  OF  SHOCK  TUNNEL  NOZZLE  THROAT  AREA 

In  order  to  size  a  suitable  nozzle  to  be  driven  by  the 

high  pressure  and  temperature  in  the  reflected  region,  a 

simplified  model  was  developed.   The  model  was  based  on  the 

following  assumptions: 

A  constant  volume  reservoir  corresponding  to  the  in- 
itial reflected  shock  conditions  was  available  to  drive 
the  nozzle  flow.   The  reservoir  volume  was  defined  by 
the  diameter  of  the  tube  and  the  axial  distance  from 
the  endplate  of  the  tube  to  the  point  where  the  contact 
surface  and  reflected  shock  met  (Fig.  3); 

The  gas  remaining  in  the  reflected  region  expands 
isentropically  to  fill  the  constant  volume; 

The  flow  through  the  nozzle  throat  was  isentropic  and 
always  choked.  (This  implies  that  the  evacuated  dump 
chamber  had  sufficient  volume,  so  as  not  to  influence 
the  nozzle  flow.) 

The  total  mass  initially  in  the  reflected  reservoir  is 

given  by, 


mo  ■  rh  vo  -  <B-1) 

g  5 


where 


2 

ttD  x 


The  mass  at  any  later  time  is  given  by 
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m  = 


R  T. 
g  t 


V, 


(B.2) 


As  the  mass  flows  through  the  choked  nozzle  the  rate  of 
change  of  mass  in  the  volume  is  equal  to  that  which  flows 
out  through  the  nozzle,  thus 


dm 
dt 


A*P, 
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Integrating 


t         t  A*P 

/   dm  =   / C,  dt 

0         0     /t7 


(B.4 


Since  the  expansion  in  the  reservior  is  assumed  to  be  isent- 
ropic,  the  stagnation  pressure  and  temperature  are  related 
according  to 
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Taking  the  ratio  of  Equations  B.l  and  B.2  and  substituting 
into  Equation  B.5, 
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Substitution  into  Equation  B.3  gives 


dm 
dt 
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so  that 
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On  integrating,  the  variation  of  mass  in  the  reflected 
reservior  with  time  is  given  by 


m0  ^  /   2  \ 
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and  the  variation  of  pressure  in  the  reflected  region  with 
time  is  given  by 
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To  estimate  a  suitable  nozzle  throat,  Equation  B.10  was 
programmed  and  a  ratio  of  Pt/P5  equal  to  one  half  was 
selected  as  the  criterion  for  calculating  the  useful  run 
time.   Using  throat  area  as  a  variable  all  other  parameters 
were  a  function  only  of  the  incident  shock  strength  and  of 
the  physical  size  of  the  tube.   Thus  Equation  B.10  enabled 
the  available  time  for  steady  flow  through  a  given  throat 
diameter  to  be  estimated. 


83 


APPENDIX  C 
ESTIMATION  OF  THE  REQUIRED  DUMP  CHAMBER  VOLUME 

In  a  shock  tunnel  application,  the  dump  chamber  pressure 
is  required  to  remain  below  the  nozzle  exit  static  pressure 
for  the  useful  duration  of  the  ref lected-shock  supply  condi- 
tions in  the  shock-tube.   To  aid  in  sizing  a  dump  chamber 
for  a  Naval  Postgraduate  School  shock  tunnel  the  following 
analysis,  based  on  a  simplified  model,  was  developed.   The 
model  assumes  the  pressure  and  temperature  are  uniform  in 
the  dump  chamber,  and  that  the  velocity  there  is  effectively 
zero.   It  is  assumed  that  the  flow  to  the  exit  of  the  nozzle 
is  isentropic,  and  that  the  volume  of  the  nozzle  is  negligi- 
ble in  comparison  to  that  of  the  chamber.   The  volume  of  the 
dump  chamber  is  given  in  terms  of  the  initial  mass  contained 
and  its  properties,  by 


RT  m 

y^  =   D  D  (C.l) 

VD      p 
D 


Dividing  by  the  reflected  reservoir  volume,  similarly  writ- 
ten in  terms  of  the  initial  mass  contained  and  its  proper- 
ties, the  ratio  of  dump  chamber  to  reflected  shock  volume  is 
given  by 


v^ 

T^ 

Pi-       HU 

D 

D 

5       D 

vrt 

Tc 

P„   m„ 

0 

5 

D       0 

(C.2) 
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As  the  flow  continues,  the  original  reflected  pressure  (P5) 
decreases  to  lower  values  (Pt)  .   Since  the  nozzle  area  ratio 
is  fixed,  the  pressure  at  the  exit  of  the  nozzle  (Pe)  also 
steadily  decreases.   Meanwhile,  the  pressure  in  the  dump 
chamber  (PD)  steadily  increases  until,  at  a  point  in  time, 
PD  =  Pe,  and  the  useful  test  time  is  over.   The  ratio  of  the 
initial  reflected  pressure  to  dump  chamber  pressure  is  then 
given  by 

_J>  =   _£  •  -±   =  (-1)   (_E)  (C.3) 

p    p     P     VP    P 
D     t     D      t    e 


where 


Y 


-t  =  1  +  -V-  M  2  (C-4) 

p  2    e 

e 


A  mass  balance  gives 


mD  =  m0  *  mi  +  mD0 


so  that 


Using  Equation  B.6  from  Appendix  B 


(C.5) 


m^       m.    m~n 

_£  =  i  -  -±  +  -2°  (c.6) 

mo      mo   mo 
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m. 

1 

m0 


P   1A5 


(C.7) 


and  substituting  in  Equation  C.6, 


^D         .        /M1/Y5+^D0 

m0  VP5/  m0 


(C.8) 


Using  Equation  C.2  applied  to  the  initial  conditions, 


m 


DO 


m, 


VUpdc 


DO 


(C.9) 


Combining   Equations    C.2    through   C.9 


3 


D 


1  + 


fe> 


Tr 


LD0 


D 


1/Yt 


D 


Yc-1 


1  + 


M 


^S"1 


(CIO) 


Equation  CIO  provides  a  relationship  between  VD/V0  and  the 
value  of  Pt/P5  which  terminates  the  useful  test  time,  in 
terms  of  given  initial  reflected  reservoir  and  dump  chamber 
conditions.   Using  an  interval  halving  technique  the  dump 
chamber  volume  requirements  were  investigated  using  an 
assumption  that  the  temperature  in  the  dump  chamber  was 
equal  to  its  initial  value,  and  also  with  the  worst-case 
assumption  that  it  was  equal  to  T5.   (It  is  noted  that  the 
temperature  change  can  be  included  explicitly  through  the 
conservation  of  energy.   However  the  complexity  of  doing  so 
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is  not  warranted  when  an  estimate  of  the  needed  volume  is 
all  that  is  required.) 
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APPENDIX  D 
SUMMARY  OF  EXPERIMENTAL  DATA 

This  appendix  contains  data  taken  using  thin  film  gauges 
and  the  pressure  transducer.   First,  in  Figure  42,  the 
measurements  obtained  of  shock  Mach  number  using  oscillo- 
scopes (preceeding  Run  50)  are  shown.   For  all  runs  from  Run 
50  on,  data  were  taken  using  the  Scientific  Atlanta  SD380 
signal  analyzer.   Tables  6  and  7  list  data  for  runs  using  a 
mixture  of  helium  and  nitrogen  and  for  runs  using  pure 
helium  as  a  driver,  respectively.  For  each  run,  trace  A 
recorded  the  response  of  the  thin  film  at  location  A,  trace 
B  recorded  the  response  of  the  thin  film  response  at  station 
B,  trace  C  recorded  the  response  of  the  endplate  pressure 
transducer.   Delay  times  were  programmed  prior  to  each  run. 
Traces  B  and  C  were  delayed  by  the  same  time.   Figure  43 
provides  definitions  of  items  which  appear  on  the  SD380 
screen  and  identifies  the  menu  selections  of  interest. 

The  last  35  pages  of  the  appendix  contain  printouts  of 
the  SD380  screen  for  the  runs  listed  in  Tables  6  and  7. 
Each  printout  can  be  interpreted  with  reference  to  Figure 
43.   Times  were  obtained  for  shock  travel  between  trans- 
ducers on  traces  A,  B  and  C,  using  the  time  scale  given  on 
the  printout,  and  the  delay  times  listed  in  Tables  6  and  7. 
These  results  are  included  in  the  two  tables. 
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Oh 


LEGEND 

n     IIK/Nli  1KEAI. 

o     MEASURE  TIM]']  A  TO  ft 
a     MEASURED  TIME  A  TO  0 
i     MEASURER  TIME  II  TO  C 


OSCILLOSCOPE  data 


MAOII  // 


Figure  42. 


Shock  Strength  vs.  P./p,  for  Data  Reduced 
from  Polaroid 
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TABLE  6 

TIME  MEASUREMENT 
(70%  HELIUM  30%  NITROGEN  DRIVER/N2  DRIVEN) 


Tempera- 

Delay 

Delay 

Time  A 

Time  B 

Time  C 

Run 

ture 

A 

B 

to  B 

to  C 

to  B 

# 

(F) 

(msec) 

(msec) 

(msec) 

(msec) 

(msec) 

62 

63 

-1.01 

0.00 

2.84 

1.67 

3.63 

63 

65 

-1.01 

-1.01 

1.77 

1.77 

1.72 

69 

64 

-1.01 

0.49 

1.81 

— 

— 

70 

64 

-1.01 

0.51 

1.36 

1.83 

3.86 

71 

64 

-1.01 

0.51 

1.72 

2.03 

3.79 

72 

65 

-1.01 

0.51 

1.76 

2.03 

3.82 

73 

64 

-1.01 

0.51 

1.80 

2.07 

3.86 

TABLE  7 

TIME  MEASUREMENT 
(HE  DRIVER/N2  DRIVEN) 


Tempera- 

Delay 

Delay 

Time  A 

Time  B 

Time  C 

Run 

ture 

A 

B 

to  B 

to  C 

to  B 

# 

(F) 

(msec) 

(msec) 

(msec) 

(msec) 

(msec) 

50 

66 

0.00 

0.49 

1.45 

— 

51 

66 

0.00 

0.49 

1.64 

0.94 

2.50 

52 

62 

-1.01 

0.49 

0.92 

1.24 

2.48 

53 

62 

-1.01 

0.49 

1.30 

0.92 

2.55 

57 

62 

-1.01 

0.49 

0.92 

1.09 

2.50 

58 

62 

-1.01 

0.49 

1.19 

1.13 

2.48 

59 

62 

-1.01 

0.49 

1.17 

1.23 

2.34 

60 

62 

-1.01 

0.49 

1.32 

1.11 

2.44 

61 

62 

— 

— 

1.15 

-- 

64 

64 

-1.01 

0.49 

1.23 

— 

65 

63 

-1.01 

0.49 

1.30 

— 

66 

61 

-- 

— 

— 

— 

67* 

63 

— 

— 

1.18 

— 

68* 

63 

-- 

-- 

1.18 

-- 

74 

67 

-2.02 

0.51 

0.97 

— 

75 

67 

— 

— 

— 

— 
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SETUP     GRP  TIME 
00:  16:  11   SYNC  TIME  C 


TRC  C 
IMP 


-© 


~© 


VW  80DB  CM  ABCD  FR  10KHZ 
UG  X2      WTG  H  C  5V 


'& 


(1)" 


0        NORM  L  NX  1   BASE 
0092969  SEC    C:   -.00613  V 

Menu  Selections  of  Interest 


© 

RUN  74 


AT  .039m  SEC    .020000 

TIME  SUM  N  0 


1  -  DISPLAY  TRACE  SELECTION 

2  -  FREQUENCY  RANGE 

3  -  Y  AXIS  SCALING 

4  -  DISPLAY  GAIN 

5  -  INPUT  LEVEL 


6  -  y,AXIS  UNITS 

7  -  RUN  NUMBER 

8  -  X  AXIS  UNITS 

9  -  X  AXIS  SCALING  (msec) 


Figure  43.   SD380  Menu  Selection  Items  of  Interest 
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APPENDIX  E 
COMPUTER  PROGRAM  FOR  SHOCK  TUNNEL  EVALUATION 

The  attached  program  was  used  for  the  evaluation  of  the 
NPS  shock  tube  and  shock  tunnel.   Sources  of  shock  tube 
relationships  are  noted  in  the  program  were  applicable. 
Output  was  tabulated  by  the  program  and  graphical  analysis 
was  provided  by  the  plotting  routine  'EASYPLOT'  installed  on 
the  NPS  IBM  computer. 
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FILE:  SH0CK2    FORTRAN   Al 

CC       THIS  PROGRAM  CALCULATES  THE  IDEAL  PERFORMANCE  OF  A  SHOCK 

CC       TUBE  USING  THE  ANALYSIS  PRESENTED  IN  THE  'FUNDAMENTALS 

CC       OF  GAS  DYNAMICS'  BY  J.  A.  OlICZAREK  PG  395-393  (FIG  8-22) 

C        &    'THE  SHOCK  TUBE  IN  HIGH  TEMPERATURE  CHEMICAL  PHYSICS'  BY  GRA 

C        GRAYDON  HURLE  1963. 

C 

C 

C 

C 

IMPLICIT  REAL  X  <\    (A-H,M,P-Z) 

REAL  FCH,X1,X2,XT0L,FT0L,FCN1,FCN2,FCN3,FCN4 

INTEGER  I,NLIM, J, II, 112 

EXTERNAL  FCN, FCN1 , FCN2, FCN3 , FCN4 

DIMENSION  Ml (200), M2 (200), M3 (200), P4P1(200),P2P 1(200) 

DIMENSION  P3P4(200),MS(200),US(200),T2(200),P5P2(200),P5P4(200) 

DIMENSION  C3(200),C4(200),P4P3(200),T3(200),TC(200),XC(200) 

DIMENSION  VC(200),P5P1(200),T5T1(200),UR(200),DELTIM(200) 

DIMENSION  DIST(200),TCC(200),XC1(20  0),DEL1A(20  0),PDI(8) 

DIMENSION  XMS ( 20 0),XP2P 1(200) ,XP5F2(200),XP3P4(200),XR02R1(200) 

DIMENSION  XT2T1(20  0),URUS(20  0),TSTART(200),VDVO(8),VD(8),AA(8) 

DIMENSION  P5(200),T5(200),DTHR(200),DU200),Mr(200),DTEST(200) 

DIMENSION  PT(200),TT(200),T(200,4),TMASS(200),DTDTHR(5),P1(200) 

DIMENSION  FT(4,10),PTPTO(10) 

COMMON  /B/  R1,MM,C8,DTHR,T 

COMMON  /A/  C0NST,GAMMA1,GAMMA4,A4A1SQ,C1,C1SQ,CC3,CC4,C5,A4A1,C9 

COMMON  /C/  ATASTR,MT,DTEST 

C       COMMON  /E/  DTHR,T,MT,DTHRDT,MM,C8,R1 
COMMON  /D/  P5,T5,T1,AA,PDI, J, II, 112 
COMMON  /F/  PTPTO, FT 
DATA  XI, X2, XT OL,FTOL, I, NLIM/0. 1,6 .0,0.0001,0.00001,0,50/ 

C 

0================================================================ 

c 

C       SELECT  GAS  CONSTANTS 

C 

C 

C 

C 

C 

C 

C 

C 


C 
C 
C       DRIVER  GAS  CONSTANTS 

R4=8314/28.02 

T4=288 

GAMMA4=1 .40 

P4=600*6 .8948E03 
C       A  IS  INCREMENT  ON  P4P1  LOOP 

A  =  100 
C      MM  IS  THE  NUMBER  OF  LOOPS 
MM  =  6 

P4P1(1)=20 
C  =================================================== 

c 

C       INPUT  TUNNEL  PARAMETERS( DRIVER  LENGTH  LI  &    L4) 

Ll=6.453 

L4=2.375 

DT1=3*.0254 
C      TO  EVALUATE  TEST  CONDITIONS  SET  FL2=1 

THR1=.0065 
FL1  =  0 
C       FL2  CALCULATES  TEST  SECTION  MACH  * 

FL2  =  1 
C       FL3  IS  FOR  NOZZLE  SIZING  CALCULATIONS 

FL3  =  0 
C      CALCULATES  TEST  SECTION  CONDITIONS 

FL4  =  1 
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GAS 

MM 

GAMMA 

SPEED  OF  SOUND 

AIR 

28.97 

1.4 

345  M/S 

ARGON 

39.94 

1  .67 

321 

HE 

4 

1  .67 

1062 

N2 

28.02 

1  .4 

350 

C02 

44.01 

1.29 

269 

DRIVEN 

GAS  CONSTANTS 

Rl=8314/28.02 

Tl=288 

GAMMA 1= 

=  1  .4 

Pl  =  l  .0E04 

FILE:  SH0CK2    FORTRAN   Al 

C      FL6  EXPERIMENT  OUTPUT  RUNS 

FL6  =  0 
C      FL7=1  PRINTS  MISC  RELATIONSHIPS 

FL7  =  0 
C      FL8  =  1  EVALUATES  PT/PTO  RATIOS  FOR  GIVEN  THROAT 

FL8  =  1 
C      EVALUATE  SPEED  OF  SOUND  CONSTANTS 
A1=SQRT(GAMMA1*R1XT1) 
A4=SQRT(GAMMA4*R4*T4) 
A4A1=A4/A1 
A4AlSQ=A4Alxx2 
C 
C=  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

C       EVALUATE  CONSTANTS  USED  COMMONLY  USED 
CI =( GAMMA 1+1)/(2XGAMMA1) 
C1SQ=C1**2 
C5=  (-2XGAMMA<4)/(  GAMMA4- 1 ) 
C6=(GAMMA4-1)/(2*GAMMA4) 
C7=(GAMMA4+1  )/(4*GAMMA4) 
C8=(GAMMA1-1)/(2XGAMMA1)X(1) 
C9=(GAMMA1+1)/(2X(GAMMA1-1)) 
C10  =  (GAMMA1/(GAMMA1-D) 
C= ============================================================= 

C       THIS  LOOP  INCREMENTS  PRESSURE  P4/P1  BY  A   AND  THEN 

C       USES  INTERVAL  HALVING  TO  DETERMINE  M3 . 

C 

DO  1  J=1,MM,1 
C0NST=P4P1( J) 

CALL  I!ITHV(FCN,X1,X2,XR,XT0L,FT0L,NLIM,I) 
M3(J)=XR 
C3(J)=CC3 
C4(J)=CC4 
X1  =  XR 
X2=6.0 
1  =  0 
C    THE  ROOT  HAS  BEEN  FOUND  AT  XR  IF  I  EQUALS  1  OR  2 
K  =  J  +  1 
P4P1(K)=P4P1(J)+A 
1      CONTINUE 
C 

C     THIS  LOOP  USES  P4/P1  8  M3  TO  DETERMINE  P2 , PI , MS, M2, US , P3P4 
C 

C 
C 


DO  5  1=1, MM, 1 

FIND  P2/P1  (EQN  8-57) 


P2P1(I)=P4P1(I)/C4(I) 

M2(I)=S0RT(2x(P2Pl(I)+(l/P2Pl(I))-2)/(GAMMAlx(GAMMAl-l)x 

&((GAMMA1+1)/(GAMMA1-1)+P2P1(I)) )) 
C       FIND  MS  8  SPEED  OF  SHOCK  (EQHC8-51) 

MS(I)=SQRT((P2P1(I)-1)XC1+1) 

US(I)=MS(I)*A1 
C       FIND  P3/P4 

P3P4(I)=C3(I)XXC5 

P4P3(I)=1/P3P4(I) 
C        EVALUATE  MAXIMUM  DURATION  OF  UNIFORM  FLOH  IN  THE  REGION 
C        BETWEEN  THE  SHOCK  HAVE  AND  THE  CONTACT  SURFACE. 
C       REFERENCE  PAGE  433  PROBLEM 
C 

T2(I)=(2/(P3P4(I)xxC7)-(4xA4x(l-P3P4(I)x*C6)x(l/P3P4(I))*xC7)/ 
8((GAMMA4-1)XA1XSQRT(1+(GAMMA1+1)X(P2P1(I)-1)/ 

&(2*GAMMAl))))xL4/A4 
C 
C 

C        T3  IS  THE  THEORITICAL  MAXIMUM  DURATION  OF  UNIFORM  FLOW 
C         IN  THE  REGION  BETWEEN  THE  CONTACT  SURFACE  AND  THE  TAIL  OF 
C        CENTERED  EXPANSION  WAVE  (REGION  3).  (PG  433  PBLM  8-6) 
C 

T3(I)=(L4/A4)*((GAMMA4-1)/2)X((P4P3(I)XX(3-GAMMA4)) 

&xx(l/(GAMMA4*4)))/(l-(P3P4(I)xxC6)) 
C 
C 
C         TC  &  XC  ARE  THE  TIME  AND  DISTANCE  AT  WHICH  THE  WAVE  FORMED  BY 
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FILE:  SH0CK2    FORTRAN   Al 

C         THE  REFLECTION  OF  THE  CENTERED  EXPANSION  HAVE  OVERTAKES  THE 

C         CONTACT  SURFACE.  (REF  PG433  PBLM  8-7) 

C 

TC(I)=(L4/A4)X2*(P4P3(I)**C7) 

XC(I)=(L4x4/(GAMMA4-l)X(l-(P3P<UI)XXC6)/(P4P3(I)xxC7))) 
C 

0====================================================================== 

c 

C        CALCULATION  OF  REFLECTED  REGION  PARAMETERS  ACCORDING  TO  REF  2 

C 

C  VELOCITY  OF  CONTACT  SURFACE  VC 

VC(I)=2XA1/(GAMMA1+1)X(MSCI)-1/MS(I)) 
C  PRESSURE  AND  TEMPERATURE  IN  REFLECTED  REGION  5  ( EQN  2.31  ; 

P5P1 C I )  =  C (2XGAMMA1XMSC I )xx2-(GAMMAl-l ) )/(GAMMAl  +  l  )  ) 
&*( (<3*GAMMA1-1 )xMS(I)xx2-2*(GAMHAl-l))/< (GAMMA1-1 )*MS( I )x*2+2) ) 

T5T1(I)=((2X(GAMMA1-1 )XMS( I ) x*2+( 3-GAMMA1 ) ) 
&X((3XGAMMA1-1)XMS(I)XX2-2X(GAMMA1-1))/((GAMMA1+I)xx2x[13(l)xy2)) 
C  CALCULATE  VELOCITY  OF  REFLECTED  SHOCK  UR  C EQN  2.38) 

UR(I)=US(I)x((2+(2/(GAMMAl-l))x(l/P2Pl(I)))/ 
8((GAMMA1+1)/(GAMMA1-1)-1/P2P1(I))) 
C  OBSERVATION  TIME  AND  DISTANCE  INREFLECTED  SHOCK 

DELTIM(I)=L1X(US(I)-VC(I))/(US(I ) x(UR( I )+VC( I )  )  ) 
DIST(I)=DELTIM(I)XUR(I) 
C  TIME  TO  EXPANSION  WAVE  REACH  CONTACT  SURFACE  (EQN  4.7) 

C=((GAMMA4+1)/(2X(GAMMA4-1))) 
TCC(I)=2XL4X(1+(GAMMA4-1 )*M3(  I  )/2  )  **C/A<4 
XC1(I)=TCC(I)X((2XA1X(MS(I)X¥2-1) )/ ( (GAMMA1+1 )XMS( I ) ) ) 
C  TEST  TIME  BETWEEN  CONTACT  SURFACE  AND  SHOCK 

DELTA(I)=  TCC(I)X(1-2X(MS(I)XX2-1)/((GAMMA1+1)XMS(I)XX2)) 
P5P2(I)=P5P1(I)/P2P1(I) 
5      CONTINUE 
C 

C=  =  =  =  =  =  =  =  "  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  "  =  r  =  =  =  =  =  =  =  =  =  r  =  =  r  =  =  =  =  =  =  = 

C===========================================  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

C      SECOND  METHOD  OF  CALCULATIONS 
Xl  =  .l 
X2=6.0 
XTOL=0.0001 
FTOL=0. 00001 
1  =  0 

HLIM=50 
XR=.5 

DO  2  J=1,MM,1 
C0HST=P4P1CJ) 
CALL  INTHV(FCN1,X1,X2,XR,XT0L,FT0L,NLIM,I) 
XMS(J)=XR 

XP2PK J)  =  (2XGAMMA1XXMS( J ) x*2-( GAMMA  1 -1 ))/( GAMMA  1  +  1 ) 
XP3P4(J)=1/(P4P1(J)/XP2P1(J)) 

XR02RK J)=((GAMMA1+1)XXMS(J)XX2)/((GAMMA1-1)XXMS( J)xx2+2) 
XT2TK J)=((GAMMA1XXMS(J)XX2-(GAMMA1-1)/2))X((GAMMA1-1)/2XXMS( J)XX 
82  +  1 )/(( (GAMMA1+1 )/2 )*X2*XMS( J )XX2 ) 
C  XP2Pl(J)=(2xGAMMAlxMS(J)xx2-(GAMMAl-l))/(GAMf1Al+l) 

XP5P2(J)=(((GAMMAl+l)/(GAMMAl-l))+2-l/Xr2Pl(J))/ 
8(1  +  (GAMMA1  +  1)/(GAMMA1-I)xl/XP2PKJ)) 
C  XP3P4(J)=1/(P4P1(J)/XP2P1(J)) 

C  XR02R1(J)=((GAMMA1+1)XMS( J)XX2)/((GAMMA1-1)XMS( J)xx2+2) 

C     XT2T1(J)=((GAMMA1XMS(J)XX2-(GAMMA1-1)/2))X((GAMMA1-1)/2XMS(J)XX 
C     82+l)/(((GAMMAl+l)/2)X*2xMS(J)XX2) 
X1  =  XR 
X2  =  6  .0 
1  =  0 
C       M3  HAS  BEEN  DETERMINED  BY  INTERVAL  HALVING 
C       DATA  XI, X2,XT0L,FT0L, I, NLIM/XR, 5. 0,0. 0001, 0.00001, 0,50/ 
URUS(J)=UR(J)/US(J) 
P5P4(J)=P5P1(J)/P4P1(J) 
C       P5(J)=P5P1(J)XP1 
P5(J)=P5P4(J)XP4 
P1(J)=P4/P4P1(J) 
T5(J)=T5T1(J)XT1 

WRITE(x,x)  P5P1(J),P4P1< J),P1(J),UR( J),DELTIM(J) 
2     CONTINUE 
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C  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

IF  CFL3.EQ.0)GO  TO  36 
C       THIS  SECTION  DETERMINES  P5  T5  AND  DETERMINES  TESTIME  AS  A 
C       FUNCTION  OF  NOZZLE  THROAT  DIAMETER. 
NN  =  1 
DO  3  1=1, MM 
C       REFLECTED  MO  AVAILABLE  TO  DRIVE  NOZZLE 

C=SQRT(GAMMAl/Rl)*(l-(GAMMAl-l)/2)**(-l*C9) 

TMASS(I)=P5(I)*(UR(I)*DELTIM(I))X3.1417*DT1XX2/(R1*T5(I)X4) 
C       HRITE(x,x)  P5P1CI),P5CI),T5CI),NN,TMASSCI) 

CALL  N0ZZLECP5CI),T5CI),NN,TMASSCI),C,FL3) 
C601    FORMATC 'P4P1=',D9. 3, 2X, »P1=' ,D9. 3, 2X, 'TOTAL  MASS  MO= ' , D9 . 4, 2X, 
C     &'RUN  =',I5/) 
NN=NN+1 
3      CONTINUE 

DO  309  1=1,100 

WRIT  EC  8, 6 02)  DTHRCI),TCI,1),T(I,2),TCI,3),TCI,4) 
602    FORMATC5CE9.3,2X)) 
C201    FORMAT ( » P5= ' , E9 . 2, 2X, ' T5  = ' , E9 . 2, 2X, 'NN=',I1 
C     8, 2X,  'TOTAL  MASS  M0=',E9.2) 

309    CONTINUE 
C      WRITE(8,603) 

C6  03   FORMATC/,  'THROAT  DIAM',2X,  'RUN1',6X,  'RUN2',7X,  'RUN3',7X,  'RUN4',/) 
C       WRITE(8,601)  CP4P1CI),P1CI),TMASSCI),I,I=1,4) 
C= =============================================================== 

c= =============================================================== 

36     IF  (FL8.E0.0)G0  TO  37 
C        THIS  SECTION  DETERMINES  P5  T5  AND  DETERMINES  TESTIME  AS  A 
C        FUNCTION  OF  PT/PTO  FOR  A  GIVEN  THROAT  DIAMETER 
NN=1 
DO  82  1=1, MM 
C        REFLECTED  MO  AVAILABLE  TO  DRIVE  NOZZLE 

C=SQRT(GAMMA1/R1)X(1-(GAMMA1-1)/2)XX(-1XC9) 

TMASSm=P5(I)x(UR(I)xDELTIM(I))x3.1417XDTlxx2/(RlxT5(I)x4) 
CALL  N0ZLE2(P5(I),T5(I),NN,TMASS(I),C,FL3) 
NN=NN+1 
82    CONTINUE 

DO  8  9  1=1, MM 

WRIT E( 8, 8 02)  PTPTOC I ) , FT( 1 , 1 ) , FT ( 1 , 2 ) , FT ( 1 , 3 ) , FTC  1 , 4 ) 
802    F0RMAT(5(E9.3,2X)) 
89     CONTINUE 

NRITECS ,803)  ( P4P1 ( I ) , 1=1 , MM) 
8  03    FORMATC/,  <U  ' P4P1 = ' , E9 . 3, 2X)  ) 
C= ===================================================== ========== 

37      IF  (FL2.EQ.0)GO  TO  133 
C       DETERMINE  TEST  SECTION  MACH  *  AS  FUNCTION  OF  THROAT 
C 

C  CALL  EXPANCTHR1,DT1) 

C       SUBROUTINE  EXPANC THR1 , DTI ) 
C      REAL  X1,X2,XT0L,FT0L,DT1,THR1 
C      INTEGER  I, HUM, J 

C      DIMENSION  MTC200),DTESTC200),DTDTHR(5) 

CC     COMMON  /A/  CONST, GAMMA1, GAMMA4, A4A1SQ, CI, C1SQ, CC3, CC4 , C5, A4A1 
C       COMMON  /C/  ATASTR,MT,DTEST,C9,DTDTHR 
WRITEC8,300) 

DTESTC1)=DT1 
Xl=.10 
XR  =  50 

DO  30  J=l,8 
X2=100 
XTOL=0.0001 
FTOL  =  0. 00001 
1  =  0 

NLIM=50 
ATASTR=CDTESTCJ)/THR1)XX2 
C      HRITECX,*)  THR1,DTEST(J),ATASTR,C9 

CALL  INTHVCFCN3,X1,X2,XR,XT0L,FT0L,NLIM,I) 
MTCJ)=XR 
X1  =  XR 
X2=100 
C 
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DTESTCK)=CDTESTCJ)+DT1*.5) 
DTDTHRC J)=DTESTCJ)/PT1 
PRINT  x, 'HELP' 

HRITEC 8, 301)  MT(J),DTDTHR(J),DTEST(J 
30  0    FORMATC/, 4X, 'MTESTMX, ' DTEST/DTUBE' 

8,4X,  'DTHR0AT',3X, ' DTUBE' ) 
301    F0RMATC4CF9.3),2C2X,E9.3)) 
30    CONTINUE 


),ATASTR,THR1,DT1 

,1X, 'DTEST',2X, 'ATASTR1 
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C 

C406 

C405 


^00 

i 

401 
21 


411 

412" 
460 
461 

402 


86 


22 
105 

i 
695 

i 
696 


IFCFL4 
THIS  S 
TEST  C 
OUTER 
DO  22 
HRITEC 
DO  21 
HRITEC  5( 
NRITE 
FORMATC 
FORMAT C 
CALL  TC 
HRITEC 
TSTARTC 
FORMATC 
TEST  TE 
FORMATC 
CONTIN 
DUMP  CH 
HRITE 
DO  461 
DO  460 
AA(I1)= 
AACI1)= 
PDICI1 
Xl=-10 
XR  =  50 
X2=100 
XTOL=. 
FTOL=. 
1  =  0 

NLIM=2 

HRITEC 

CALL  I 

VDVOCI 

HRITE 

X1=XR 

X2  =  10 

VDCI1 

HRITE 

FORMA 

X) 

FORMA 

CONTI 

CONTI 

HRITEC8 

FORMATC 

"TIME  F 

'TIME  T 

CONTIN 


.EQ.O) 

ECTION 

HAMBER 

LOOP  I 

J  =  1 , MM 

8,400) 

1  =  1,8 

,*)  CI 

CX,405 

2X, *MA 

6CD9.2 

ONDCMT 

8,401) 

I) 

/,  'TES 

MP'  ,3X 

9CE8.3 

UE 

AMBER 

C8,411 

112=1, 

11=1,1 

.1X11 


GO  TO  105 
EVALUATES  TEST  CONDTIO 
AREAS  GIVEN  A  FIXED  NO 

S  FOR  P4P1  CHANGE,  INNE 


NS  FOR  VARIOUS 

ZZLE  THROAT  DIAMETER 

R  LOOP  EVALUATE  CH  AREA 


)  MTCI) 
CH',6X, 
)) 

CI),P5C 
MTCI), 

T  MACH' 
, 'MS ',4 
,1X),E6 

SIZING 

) 

1 


,P5CJ),T5CJ),PT( 
'P5  =',6X, 'PT',5 

J),T5CJ),PTCI),T 
P5CJ),PTCI),T5CJ 

,5X, 'P5',4X, 'TES 
X,  'P4P1',6X,  'PI' 
.1) 

EVALUATION 


I),TTCI),C10 

X,'T5',8X,'TT',5X,'TSTART') 

rCI),TSTARTCI),THRl,C10) 
),TTCI),MSCJ),P4P1CJ),P1CJ) 

T  P'  ,4X, 'T5' ,4X, 
,7X, 'TSTART') 


)=. 013332 
000 

000 

0001 

000001 

00 

8,*)  AACI1),P 

NTHV(FCN4,X1, 

1  )=XR 

C8,x)  AACI1),PDICI1),VDV0CI1) 


DICI1) 
X2,XR,XTOL,FTOL,NLIM, I) 


)=VDV0(I1)*CU 
(8,412)  P5CJ) 
TC/,3X, 'P5',8 

TC5CE9.2,1X)) 

NUE 

HUE 

,402)  DELTIMC 

/,  'TEST  TIME 

OR  .5*P5  OUT 

0  START  NOZZL 

UE 


R(J)*DELTIMCJ)*3 
,T5C J),VDV0CI1), 
X,  'T5',8X, 'VDVO' 


J), TCI, J) 
OF  SHOCK  TUBE  = 
OF  NOZZLE  =',D9. 
E  CEMPERICAL  EST 


.1417*DTl**2/4) 

VDCI1 ),AACI1) 

,7X, »VD',8X, 'PTPTO', 


',D9.4,/ 
4) 
IMATE)',D9.4) 


K=MM 

IF    CFL6.EO.0)GO    TO   699 
HRITEC8,695) 

HRITEC 8, 6 96)    C P4P1 ( J ) , MSC J ) , USC J 
,P5CJ), J=1,K) 

FORMATC 2X, 'XP4P1',7X, 'MS',9X,  'VS' 
5X, 'P5*,8X, 'T5') 

F0RMATC6CF9.3,2X),E9.3) 


),VCCJ),URUSCJ),T5CJ) 
,9X,  'VC9X,  'URUS', 


699         IF    (FL1 .EQ.0)GO    TO    108 

HRITEC8,201)  131 
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NRITE(8,200)  ( PAP1 ( J ) , MS( J ) , XP2P1 ( J ) , XP5P2( J ) , XP3PA( J ) 
8,XR02R1(J),XT2T1(J),XMS(J),URUS(J), J=1,K) 

200  F0RMAT(9(F11.3,2X)) 

201  FORMAT (2X,  'XPAP1',7X,  'MS!,8X,  'XP2P1»,8X,  'XP5P2',6X,  'XP3PA 
&,AX, 'R0U2/R0U1',AX, 'XT2/T1',6X, 'XMS',6X, 'UR/US' ) 

MRITEC8, 203)  ( PAP1 ( J ) , MS( J ) , XT2T1 ( J ) , T 5T1 ( J ) , P2PK J ) 
&,P5PA(J),DELTIM(J),J=1,K) 
203      FORMAT(7(D9.A,2X)) 
108       MRITE(8,20A) 

l!RITE(8,205)  R1,T1,GAMMA1,RA,TA,GAMMAA 
20  A      FORMAT  (3X,  ,PAP1,,8X,,MSI,8X,,XT2T1,,7X,,T5T1,,7X,,P2P1' 

&,7X, 'P5PA',5X, 'DELTIME') 
205    FORMAT ( ' Rl  =  • , F9 . 2, 2X, 'Tl  =  ' , F9 . 2,  2X, ' GAMMA  1  =  ' , F9 . 2, 
8/,  'RA=',F9.2,2X,  'TA=',F9.2,2X,  •  GAMMAS  =  '  ,  F9  .  2  ) 


C======== ============================= ======= 

IF  (FL7.EQ.0)  GOTO  109 
C        OUTPUT  DATA  TO  A  FILE 
NRITE(8,103) 

NRITE(8,100)  (PAP1(I),MS(I),M2(I),M3(I),AAA1, 
&P5P1(I),T5T1(I),P3PA(I),P2P1(I),P5P2(I),DELTIM(I),I=1,K) 
C        HRITE(8,106) 

C        WRITE(8,107)  (PAP1(I),UR(I),DIST(I),XC1(I),TCC(I),DELTIM(I), 
C     &DELTA(I),T2(I),I=1,K) 

C    106         FORMAT(AX, 'PAP1',2X, 'VEL    REFL  '  ,  AX,  •  DIS    REFL',AX, 
C  S'EXP    POS',AX,'EXP    TIME', AX, 'DEL    TIM ', 5X, ' DELTA ', 8X, ' T2 ' ) 

C  HRITE(8,10A) 

C    107  FORMAT(8(F9.A,2X)) 

C  WRITE(8,105)    (  P^iPl  ( I ) ,  T2(  I ) ,  T3(  I) ,  TC(  I  ) ,  XC(  I ) ,  1  =  1 ,  K) 

100         F0RMAT(11(F9.A,1X)) 
C    102         FORMAKAX,  'PAPl'^X,  'T2',8X,  'PSPA'^X,  'P2P1') 
103         FORMAT (AX, '  PAP1',6X, 'MS', 

89 X,  'M2' ,6X,  M13',8X,  'AAA1' ,6X,  'P5P1' 
&,6X, 'T5T1',6X, 

&1X,'P3PA',6X,,P2P1,,6X,IP5P2,,6X,'TESTTIME') 
C  10A    F0RMAK2X,  »PAP1',12X,  'T2',8X,  'T3M0X,  'TC',9X,  'XO 
C  105    FORMAT(5(F10.A,2X)) 
109    X=l 
STOP 
END 
C 
C 
C 
C 
C 

REAL  FUNCTION  FCN(X) 
REAL  X 

COMMON/ A/  CONST, GAMMA  1, GAMMA A, AAA1SQ, CI, CI SQ,CC3, CCA, C5,AAA1,C9 
CC3=(l+(GAMMAA-l)XX/2) 

CCA=(CC3*X(2*GAMMAA) )**(1/(GAMMAA-1)  ) 
FCN=(l+( . 5*AAA1SQ)*( ( GAMMA 1XX/CC3 ) **2 ) 
&*(C1+SQRT((C1SQ+A/(AAA1SQX(GAMMA1*X/CC3)*X2))))) 
JXCCA-CONST 
C     WRITE(8,200)  FCN, GAMMA1 , GAMMAA , AAA1 SQ, CC3, CCA , CONST 
C  200    F0RMAT(7(F8.2)) 
RETURN 
END 
C 

REAL    FUNCTION    FCNKY) 
REAL    Y 
COMMON/ A/    C0NST,GAMMA1,GAMMAA,AAA1SQ,C1,C1SQ,CC3,CCA,C5,AAA1,C9 
FCN1=((2*GAMMA1*Y**2-GAMMA1+1)/(GAMMA1+1))*(1-((GAMMAA-1)/(AAA1 
&X(GAMMA1+1))X(Y-1/Y) ))X*C5-C0NST 
C  PRINT*, FCN1, CONST, AAA1,C5 

RETURN 
END 


REAL    FUNCTION    FCN3(Z) 
REAL    Z 

COMMON   /C/    ATASTR,MT,DTEST 

COMMON   /A/    CONST, GAMMA1, GAMMAA, AAA1SQ, CI, CI SQ,CC3, CCA, C5,AAA1,C9 
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C     WRITE(x,x)  ATASTR,GAMMA1,Z,C9 

FCN3=ATASTR-1/ZX(2X(1+(GAMMA1-1)/2XZXX2) 
8/((GAMMAl+l)))**C9 
C     WRITE(X,x)  ATASTR,GAMMA1,Z,FCN3 
RETURN 
END 
C 

REAL  FUNCTION  FCN4CM) 
REAL  W 
C      INTEGER  J, 12, II 

DIMENSION  P5(200),T5(200),AA(8),PDI(8),MT(200) 

COMMON/C/  ATASTR,MT,DTEST 
COMMON/ D/  P5,T5,T1,AA,PDI, J,I1,II2 

FCN4=W-(l+PDI(Il)*T5(J)*W/(P5(J)xTl)-AA(Il)xx.714)XTlx 
&(1+.2*100)**3.5/(AA(I1)*T5(J)) 
C      WRITE(8,131)  AA(I1),PDI(I1),P5( J),W,T5(J), J,I1 
C  131  FORMATC 'A=',E9.2, 'PDI=',E9.2, *P5=',E9.2, 
C     &'W=',E9.2,1X,E9.2,1X,F<+.1,1X,F4.1) 
RETURN 
END 
C 
C 

SUBROUTINE  NOZZL E( X, Y, NN, MTOTAL , C0N3T2, FL 3 ) 
C      THIS  ROUTINE  DETERMINES  TIME  FOR  NOZZLE  FLOW  FOR  A  GIVEN  THROAT 
C       DIAMETER. 

REAL  X8  X,Y, MTOTAL, C0NST2, PI 
INTEGER  NN,FL3 

DIMENSION  DTHR(200),T(200,4) 

COMMON  /A/  C0NST,GAMMA1,GAMMA4,A4A1SQ,C1,C1SQ,CC3,CC4,C5,A4A1,C9 
C       COMMON  /B/  P5,T5,R1,MM,C8,DTHR,DTHRDT,T 
COMMON  /B/  R1,MM,C8,DTHR,T 
PI=3.1416 
DTHR(1)=.0065 
DO  30  1=1,100,1 

T(I,NN)=MT0TALXDSQRT(Y)X2X(1/.5XXC8-1  )/ ( C0NST2XPlxDTHR( I )**2/4*X 
&X(GAMMA1-1)) 
J  =  I  +  1 

DTHR(J)=DTHR(I)+.0005 
C      URITEC8, 20  0 )  DTHRC I ) , T ( I , NN ) 
C200    F0RMAT(2(E9.3,2X)) 

C201    FORMAT ( • P5= ' , E9 . 2, 2X, ' T5= * , E9 . 2, 2X, 'NN=',I1 
C     8, 2X,  'TOTAL  MASS  MO=',E9.2) 

30     CONTINUE 
C      MRITE(8,201)  X, Y, NN, MTOTAL 
RETURN 
END 
C 
C 

SUBROUTINE  NOZL E2( X, Y, NN, MTOTAL , CONST 2, FL 3  ) 
C       THIS  ROUTINE  DETERMINES  TIME  FOR  NOZZLE  FLOW  FOR  A  GIVEN  THROAT 
C       DIAMETER  WITH  VARATION  OF  PT/PTO 
REAL  X8  X,Y, MTOTAL, C0NST2, PI 
INTEGER  NN,FL3 

DIMENSION  DTHR(200),T(200,4),PTPTO(10),FT(4,10) 

COMMON  /A/  CONST, GAMMA  1, GAMMA4, A  4 A 1SQ, CI, C1SQ, CC3, CC4, C5, A4 A 1,C9 
C      COMMON  /B/  P5,T5,R1,MM,C8,DTHR,DTHRDT,T 
COMMON  /B/  R1,MM,C8,DTHR,T 
COMMON  /F/  PTPTCFT 
PI=3.1«U6 
DTHR(1)=.0065 
PTPT0(1)=.2 
DO  80  1  =  1,  <\ 

FT(I,NN)=MT0TALXDSQRT(Y)X2X(1/PTPT0(I)XXC8-1) 
&/(C0NST2XPIXDTHR(l)XX2/4XX 
SX(GAMMAl-l)) 
J  =  I  +  1 

PTPT0(J)=PTPT0(I)+.2 
80     CONTINUE 
RETURN 
END 
C 
C      SUBROUTINE  EXPANC THR1 , DTI ) 
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C  REAL  X1,X2,XT0L,FT0L,DT1,THR1 

C  INTEGER  I , NL IM, J 

C  DIMENSION  MT(200),DTEST(200),DTDTHR(5) 

C  COMMON  /A/  CONST, GAMMA1, GAMMA4, A4A1SQ, CI, C1SQ, CC3, CC4, C5, A4A1 , C9 

C  COMMON  /C/  ATASTR,MT,DTEST,C9,DTDTHR 

C  DTEST(1)=DT1 

C  Xl=1.0 

C  XR  =  2 

C  DO  30  J=l,2 

C  X2=10 

C  XTOL=0.001 

C  FTOL=0.0001 

CC  1  =  0 

C  HLIM=50 

C  ATASTR  =  (DTEST( J)/THR1  )**2 

C  NRITEC*,*)  THR1,DTEST( J),ATA5TR 

C  CALL  INTHV(FCN3,X1,X2,XR,XT0L,FT0L,NLIM,I) 

C  MT(I)=XR 

C  X1=XR 

C  K=J+1 

C  DTEST(K)=(DTEST(J)+DT1) 

C  DTDTHR(J)=DTEST(J)/DT1 

C  WRITE(8,301)  MT(J),DTDTHR(J),DTEST(J) 

C30  0    FORMAT C2X, 'MTEST',2X, ' DTEST/DTUBE ' , 2X, 'DTEST') 

C  1JRITE(8,300) 

C301  F0RMAT(3(F9.2)) 

C  30    CONTINUE 

C  RETURN 

C  END 

C 

C 

c 

SUBROUTINE  TCONDC Zl , P5, T5 , PT , TT , TIME, THR1 , CI  0 ) 

REAL  *8  P5,T5,TT, TIME, Z1,THR1,B, CIO 

COMMON  /A/  CONST, GAMMA  1 , GAMMA4 , A4A1SQ, CI, CI SQ, CC3, CC4, C5, A4A1 , C9 

C      COMMON  /E/  DTRH(200),T(200),MT,DTHRDT,MM,C8 

COMMON  /B/  R1,MM,C8,DTHR,T 

WRITE(X,403) 

HRITE(*,404)  Z1,P5,PT,T5,TT,C10 

HRITE(*,x)  CIO 

A5=SQRT(GAMMA1*R1XT5) 

PT=P5/(1+(GAMMA1-1)/2*Z1XX2)XXC10 

TT=T5/(1+(GAMMA1-1)/2XZ1XX2) 

C        B  IS  THE  ANGLE  OUT  OF  THE  THROAT  30  DEGREES 

B=TAN(15*3. 1416/180) 

C  TIME=0 

TIME=.76XTHR1/(2XBXA5)X((DSQRT(Z1))XQ-1/Z1) 

8+l/15x((DSQRT(Zl))XX3-l)+l/25X((DSQRT(Zl))*X5-l)) 

4  03    FORMAT (2X,  'MACH' ,6X,'P5  =I,6X,,PTI,5X,'T5I,8X,,TT,,5X,  'TSTART'  ) 

HRITE( *,4  04)  Z1,P5,PT,T5,TT, CIO, TIME 

404    F0RMAT(7(D9.2)) 

RETURN 

END 
c    

C 
C 

c 
c 
c 

SUBROUTINE  INTHV( FCN, XI , X2, XR, XTOL , FTOL , NL IM, I ) 
C 
c    

C       THIS  ROUTINE  MAS  DEVELOPED  BY  PROF.  MILLER  DEPT  AERONAUTICS 

C       SUBROUTINE  INTHV  i 

C  THIS  SUBROUTINE  FINDS  THE  ROOTS  OF  A  FUNCTION  , 

C  F(X)  =  0,  BY  INTERVAL  HALVING. 

C 

c  

C 

C  PARAMETERS  ARE  : 

C  FCN        -FUNCTION  THAT  COMPUTES  VALUES  FOR  F(X).MUST  BE  DECLARED 

C  EXTERNAL  IN  CALLING  PROGRAM. IT  HAS  ONE  ARGUMENT,  X. 
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XI, X2  -INITIALS  VALUES  OF  X.F(X)  MUST  CHANGE  SIGN  AT  THESE  PTS 

XR  -RETURNS  THE  ROOT  TO  THE  MAIN  PROGRAM. 

XTOL,FTOL  -TOLERANCE  VALUES  FOR  X,F(X)  TO  TERMINATE  ITERATIONS. 

NLIM  -LIMIT  TO  NUMBER  OF  ITERATIONS. 

I  -A  SIGNAL  FOR  HOM  ROUTINE  TERMINATED. 

1=1  MEETS  TOLERANCE  FOR  X  VALUES. 

1=2  MEETS  TOLERANCE  FORF(X) 

I=-l  NLIM  EXCEEDED. 

I=-2  F(X1)  NOT  OPPOSITE  IN  SIGN  TO  F(X2) 
WHEN  THE  SUBROUTINE  IS  CALLED, THE  VALUE  OF  I  INDICATES  WHETHER  TO 

PRINT  EACH  VALUE  OR  NOT. 1=0  MEANS  PRINT  THEM, I. HE. 0  MEANS  DON'T. 


REAL  FCN,X1,X2,XR,XT0L,FT0L 
INTEGER  NLIM, I, J 
REAL  F1,F2,FR,XERR 


CHECK  THAT  F(X1)  8  F(X2)  DIFFER  IN  SIGN 


Fl  =  FCN(Xl) 
F2=  FCNCX2) 
IF  (F1XF2  .GT 

I  =  -2 

PRINT  201 

RETURN 
END  IF 


0.0)  THEN 


COMPUTE  SEQUENCE  OF  POINTS  CONVERGING  TO  THE  ROOT 


DO  20  J=1,NLIM 

XR  =  (Xl+X2)/2.0 

FR  =  FCN(XR) 

XERR  =  ABSCX1-X2V2.0 

IF  (I  .EQ.  0.0)  THEN 

PRINT  199,  J,XR,FR 
END  IF 


CHECK  ON  STOPPING  CRITERIA 

IF  (  XERR  .LE.  XTOL)  THEN 

1  =  1 

PRINT  202,  J,XR,FR 

RETURN 
END  IF 


IF  (  ABS(FR)  .LE.  FTOL  )  THEN 
1  =  2 

PRINT  203,  J,XR,FR 
RETURN 
END  IF 


COMPUTE  NEXT  X 

IF  (  FR*F1 
X1  =  XR 
F1  =  FR 
ELSE 
X2  =  XR 
F2  =  FR 
END  IF 
20       CONTINUE 


GT.  0.0  )  THEN 


135 


FILE:  SH0CK2    FORTRAN   Al 


c    

C 

C    WHEN  LOOP  IS  NORMALLY  COMPLETED, NL IM  IS  EXCEEDED 

C 

I=-l 

PRINT  200,  NLIM,XR,FR 

RETURN 


199  FORMATC  AT  ITERATION1  ,  13,  3X,  ■  X=  '  ,  E12  .  5,  4X,  ■  F(X)  =',E12.5) 

200  FORMATC/'  TOLERANCE  NOT  MET  AFTER  ',14,*  ITERATIONS   X  =  ', 
1  E12.5, •  F(X)  =  ',E12.5) 

201  FORMATC/'  FUNCTION  HAS  THE  SAME  SIGN  AT  X1&X2') 

202  FORMATC/'  TOLERANCE  MET  IN  '  , 1 4 ,  '  ITERATIONS   X  =',E12.5, 
1  '  FCX)  =  »,E12.5) 

203  FORMATC/'  F  TOLERANCE  MET  IN   ',14,'  ITERATIONS  X  =  »,E12.5, 
1  '  FCX)  =',E12.5) 

END 

C 

c    
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